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ABSTRACT 

We present properties of individual and composite rest-UV spectra of continuum- and narrowband- 
selected star-forming galaxies (SFGs) at a redshift of 2 < z < 3.5 discovered by the MUSYC collab- 
oration in the Extended Chandra Deep Field-South. Among our sample of 81 UV-bright SFGs, 59 
have R<25.5, of which 32 have rest-frame equivalent widths WL yo > 20A, the canonical limit to be 
classified as a Lya-emitting galaxy (LAE). We divide our dataset into subsamples based on proper- 
ties we are able to measure for each individual galaxy: Lya equivalent width, rest-frame UV colors, 
and redshift. Among our subsample of galaxies with R<25.5, those with rest-frame Wl Y q > 20A 
have bluer UV continua, weaker low-ionization interstellar absorption lines, weaker C IV absorption, 
and stronger Si II* nebular emission than those with Wl yq < 20A. We measure a velocity offset of 
Av ~ 600 km s _1 between Lya emission and low-ionization absorption, which does not vary sub- 
stantially among any of our subsamples. We find that the interstellar component, as opposed to the 
stellar component, dominates the high-ionization absorption line profiles. We find the low- and high- 
ionization Si ionization states have similar kinematic properties, yet the low-ionization absorption is 
correlated with Lya emission and the high-ionization absorption is not. These trends are consistent 
with outflowi ng neutral gas being in the form of neutral clouds embedded in ionized gas as previously 
suggested by iSteidel et al.l (|2010f ). Moreover, our galaxies with bluer UV colors have stronger Lya 
emission, weaker low-ionization absorption and more prominent nebular emission line profiles. From 
a redshift of 2.7 < z < 3.5 to 2.0 < z < 2.7, our subsample of galaxies with WL ya < 20 A show no 
significant evolution in their physical properties or the nature of their outflows. Among our dataset, 
UV-bright galaxies with W Lya > 20A exhibit weaker Lya emission at lower redshifts, although we 
caution that this could be caused by spectroscopic confirmation of low Lya equivalent width galaxies 
being harder at z^3 than z^2. 



1. INTRODUCTION 

Narrowband and broadband selection techniques 
around the Lya emission line and Lyman limit respec- 
tively have made isolating Lya emitters (LAEs) and 
UV-bright, continuum-selected galaxies quite effective 
at redshifts up to z — 7. Through stacking photometric 
samples and stellar population modeling, numerous 
groups have uncovered a wealth of information on 
these galaxies' properties including stellar masses, star 
formation rates, ages, amount o f dust (E(B — V)) 
and star formation histories (e.g., [A dclbcrgc r fc Steide 
20001: iPapqvich et all I200H IShaplev et al.l 120011 ' " 
Erb et al. I2006bl iReddv et all 
2006ri IPirzkal et al l 120071 



2001 iNilsson et al.1 [20091 iFinkelstein etafl [200l 
Guaita et al.ll2010h . As a result, there have been numer- 



ous follow-up r est-frame ultraviolet ( UV) spect r oscopic 
surveys (e.g . . IShaplev et al.l 



iKornei et all 120101: iTalia et al.l 



2003 



2011 



20061: Gawiser et al 



iGawiser et all 12007 



iteidell 
2005; 
et al.1 
l2007t 



Nilsson et al.l 120071: iPentericci et al.l 120071: lLai et al l 
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lErb et al.l 12006b: 
Uones et aLll201lD . 

Rest-frame UV spectroscopy offers the ability to not 
only confirm redshifts, but also to study galactic-scale 
outflows, stellar photospheres and dust reddening 
through interstellar absorption line kinematics, the 
shape and strength of the Lya emission line, P Cygni 
profiles from ma ssive stars, and th e UV spectral slope 
ilLeitherer et al.l 119951: fPettini et al.l 120001: IShaplev et al.l 
2003). 

In order to determine the impact of feedback from star 
formation on galaxy evolution, we have targeted the peak 
of cosmic star formation density occuring at a redshift of 
1.5 < z < 3 when galaxies and their surrounding gas 
were evolving rapidly. Understanding the relation be- 
tween star-forming galaxies and the surrounding inter- 
galactic medium (IGM) at these redshifts will give insight 
into this relation at higher redshifts where spectroscopy 
is much more difficult. It is well known that outflows 
can occur in starburst galaxies in the local universe and, 
more recently, it has been shown that they a re ubiquitous 
among high-redshift star-for ming g alaxies (ISteidel et al.1 
1996at IPettini et al.l 120001 120011 : IShaplev et all 1 2003): 
Martini 120051 : iRupke et al.l l2005t iTremonti et al.l 12007 



Wcincr et al. 2009; SteioMeLalj|2010fi. Analysis of these 
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outflows has primarily been focused on modeling the 
complex spatial and velocity structures, the covering 
fracti on, and the nature o f neutral and ionized gas 
(e.g.. [Heckman et al.l 119901: iLehnert fc Heckman| 1 19961 : 
Martinlll999l: iShapley et al.ll2003t iStrickland etltlbooi 
Schw artz et al.l 2006; Grim es et al.l [2009; St eidel et al.1 



20101) . Galactic-scale outflows are not only important 



in understanding how feedback from star formation af- 
fects galaxy evolution, but also how they affect the metal 
enrichment, ionization state, and physical state of the 
IGM. 

In order to understand the connection between LAEs 
and continuum-selected galaxies, it is important to un- 
derstand the different classifications for each type of 
galaxy, namely the criteria by which they are selected. 
Continuum-selected galaxies are chosen via color cuts 
and apparent magnitude limits, while LAEs are selected 
via equivalent width and line flux cuts that are i ndepen - 
dent of their colors. As defined bv lSteidel et al.l (|T996b), 
Lyman break galaxies (LBGs) are continuum-selected 
galaxies at z~3 and typi cally have an apparen t mag- 
nitude limit of R<25.5. IGawiser et all (|2006bD found 
LAEs at z^3 with LL yQ > 4 x 10~ 17 ergs s _1 cm~ 2 to 
have a median apparent ma gnitude of R = 27. In spite 
of this luminosity difference. IGawiser et al.l (|2006bD find 
that z — 3.1 LAEs have similar rest-frame UV colors 
to LBGs, indicating that they are missed in continuum- 
selected su rveys primarily due to their fainter continua. 
Similarly, (jShaplev et al.ll2003l hereafter S03) find that 
20 — 25 percent of their almost 1000 LBGs have strong 
enough Lya emission lines (Wo > 20A) to be classified 
as LAEs. 

Given the faint nature of these galaxies (R> 23), pre- 
vious analyses have focused on stacking large numbers 
of low signal-to-noise spectra. Similar to stacking photo- 
metric samples, this technique yields the average prop- 
erties of a 'typical' galaxy. S03 applied this technique 
to almost 1000 LBGs at a redshift of z~3 and found 
that they display a wide distribution of properties, with 
Lya ranging from strong emission to absorption, varying 
UV continuum slopes, a range of velocity offsets between 
Lya emission and interstellar absorption, and different 
strengths of interstellar absorption lines. Through divid- 
ing their sample into quartiles of Lya equivalent width, 
they found that the sample with strongest Lya emission 
had bluer UV continua, weaker low-ionization absorp- 
tion, smaller kinematic offsets, stronger nebular emis- 
sion, and older stellar populations. With the exception 
of older stellar populations, these trends suggest that 
galaxies with stronger Lya emission are less evolved, less 
rapidly forming stars, and likely less massive. Nonethe- 
less, the presence of older stellar populations implies 
more evolution, indicating that galaxies with strong Lya 
emission may be a heterogeneous population. 

As the defining characteristic of LAEs, it is vital to 
understand the origin of the Lya emission line, which 
is largely dependent on the effects of radiative trans- 
fer. Lya photons are subject to resonant scattering and 
therefore have much longer pathlengths than continuum 
photons, making them more subject to extinction from 
dust. Thus, one explanation for LAEs' large Lya equiv- 
alent widths is that they are young, chemically pristine 
galaxies undergoing one of their first major episodes of 



star formation. A second hypothesis suggests that they 
are older, more evolved galaxies where dust is segregated 
into clumps that are surrounded by a relatively dust-free 
intercloud mediu m (iNeufeldl 119911 : Hansen fc Ohl 120061 : 
iFinkelstein et al.1 120091 ). In this scenario, Lya photons 
scatter off dust clumps while UV continuum photons 
penetrate them and are subject to extinction, thereby 
increasing the relative strength of the Lya emission line. 
A lack of correlation between the spectral index of the 
UV continuum slope and Lya equivalent width is inter- 
preted as evidence for the decoupling of continuum and 
line extinction (S03) as would be expected if the neu- 
tral medium was i nhomogeneous. For 2 < z < 4 LAEs, 
iBlanc et all (|2011ft recently found E(B-V) to correlate 
with Lya escape fraction in such a way that Lya photons 
are being extinguished by the same amount as UV contin- 
uum photons, leading them to suggest that among LAEs, 
the Lya equivalent width is neither being enhanced due 
to a clumpy I SM nor being preferentially quenched by 
dust. Finally, iMalhotra fc Rhoadsl (|2002l ) suggest that 
LAEs' large Lya equivalent widths may instead be due 
to a peculiar initial mass function, zero-metallicity stars, 
or narrow-line active galactic nuclei. 

Through stellar population modeling, LAEs have been 
shown to encompass a broad range of physical proper- 
ties. At an epoch of 2 < z < 4.5, many groups find that 
they are characterized as young (< 100 Myr), relatively 
dust-free (E(B-V)< 0.3) galaxies with low stellar masses 
(< 10 9 Mp)) an d low star formation rates ( 2 — 3 M^yr" 1 ) 
(IGawiser et al.ll2006bl: iPirzkal et al.ll2007t IGawiser et al.l 
l2007tlNilsson et alJl2007l ). Furthermore. LAEs have very 
high specific star formation rates (7 x 10~ 9 yr" 1 ), de- 
fined as star formation rate per unit mass, as compared 
to LBGs (3 x 10~ 9 yr -1 ) indicating they are building 
up stellar mass at a faster rate than continuum-selected 
galaxies. This picture is consistent with LAEs being 
chemically pristine and at the beginning of an evolu- 
tionary sequence. On the other hand, various authors 
have found LAEs to exhibit a broader range of proper- 
ties with ages ranging up to > 1 Gyr, stellar masses up to 
10 10 M , and A 1 200 (defined as th e absor pt ion at 1200 A) 
as high as 4.5 dPentericci et all (|2007D ; iNilsson et~aLl 
(2009,); IFinkelstein et al.l (|2009T )) These results indicate 
that some LAEs are older, dustier, more evo l ved g alax- 
ies similar to LBGs . Furthermore, lLai et afl (|2008| ) and 
iGuaita et al.l (|2010l) found IRAC detected obj ects to ex- 
hibit properties similar to LBGs, motivating lLai et afl 
(2008) to suggest they may be a lower mass extension of 
LBGs. 

In this paper, our goal is to understand the relation- 
ship among Lya emission, interstellar kinematics, mul- 
tiphase gas, nebular emission, and spectral slope for a 
sample of 2 < z < 3.5 UV-bright star-forming galaxies. 
We subdivide our sample in terms of their Lya equiva- 
lent width, rest-UV colors, and redshift to examine the 
relationships among these properties. This paper is or- 
ganized as follows. In Sj2j we present the spectroscopic 
sample including details of the observations and data re- 
duction. We discuss the spectroscopic properties of our 
individual spectra in $3] In J|] and [5j we describe the 
stacking technique and analyze the composite spectra as 
subdivided based on Lya equivalent width. Selection ef- 
fects and trends among the outflowing component and 
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galaxy properties are presented in ££SJ Finally, in fJTl we 
summarize our conclusions. 

2. OBSERVATIONS AND METHODOLOGY 

The galaxies that will be presented and analyzed in 
the following sections include narrowband-excess objects 
and continuum-selected sources. The latter category in- 
cludes U band drop-outs at z^3 and galaxies with promi- 
nent Lya decrements at z^2 (LBGs and "BX" galax- 
ies, continuum -selected UV- bright SFGs identified via 
the criteria of iSteidel etaLl I1996U I2004D . The imag- 
ing data used to select our targets were collected us- 
ing the MOSAI C II CCD camera on the CT IO Blanco 
4 m telescope (http://www.astro.yale.edu/MUSYC/). 
The narrowband data consisted ot a series ot exposures 
taken through 50A wide full width at half-maximum 
(FWHM) filters at A5007A and A4990A for z = 3.1 
LAEs and A3727A for z = 2.07 LAE s. The details 
of the se o bservations are p resen ted in IGronwall et al.l 
(20071 and ICiardullo et aU (I20H for LAEs at z = 3.1 
and IGuaita et al.l (|2010|T at z = 2.07. O ur LBG selec- 
tion criteria is presented in lGawiser et al.l (|2006al) . which 
covers a redshift interval of 2.6 < z < 3.5, wh i le our 
BX selection criteria is presented in IGuaita et al.l ()2010f ) 
and spans 2 < z < 2.7. Our spectroscopic candidates 
were chosen from these two photometric surveys. In 
this paper, we present spectra taken from two spectro- 
scopic surveys using the FORS and VIMOS instruments 
on the Very Large Telescope (VLT). The first survey, us- 
ing the FORS instrument, targeted primarily z^2 LAEs 
and BXs, along with a few z~3 LAEs and LBGs. The 
primary targets for observation using t he VIMOS in- 
strum ent were X-ray sources analyzed in iTreister et al.1 
(2009) with z-3 LAE and LBG candidates used to fill 
the remaining slitlets. 

2.1. FORS Spectroscopy 

For our FORS spectroscopic sample, the primary tar- 
gets were chosen by first prioritizing z = 2.1 LAE candi- 
dates, followed by BX candidates and finally z^3 LAEs 
and LBGs were used to fill available slitlets. The targets 
covered six 7' x 7' masks in the CDF-S, however one of 
these masks was unusable due to poor observing condi- 
tions and has been omitted from the dataset. Each of the 
five remaining slitmasks consisted of ~ 30 targets with 
~ 70% expected to be at z^2. 

These data were obtained using the FORS2 spectro- 
graph on UT1 with the B1200 grism, the c2V UV- 
scnsitive CCD, an atmospheric dispersion corrector, and 
l."0 slitlets. This configuration led to a dispersion of 
0.86A per pixel corresponding to a nominal spectral res- 
olution of 3.4A (FWHM) or R~1400 over a spectral cov- 
erage of 3200 — 6000A for each slitlet. The small angular 
size of our targets caused the actual spectral resolution 
to be seeing limited in almost all cases, which varied from 
0."7 to l."2 leading to a minimum spectral resolution of 
2.4A. The observations were taken in a series of four ex- 
posures with exposure times of 1800 — 2400 seconds per 
mask. We also note that the throughput at A3727A is 
~2 times worse than that at A5000A. 

The data were reduced using the standard FORS 
pipeline. This process included bias-subtraction, flat 



fielding of science data, sky subtraction, and wave- 
length calibration through arc lamp exposures. For 
5 objects, we had difficulty wavelength calibrating the 
A3200 - 4200A region. We found an optimal wavelength 
calibration by lowering the threshold for identifying arc 
lamp lines and adding two lines to the arc lamp file at 
3612.87A (Cd I) and 3663.00A (Hg I). We removed three 
objects that had data in this region for which we were 
not able to correctly calibrate their wavelengths. The 
spectra were then divided into individual spectrograms 
whereby one dimensional spectra were extracted using 
IRAF's apall task. Finally, each one-dimensional spec- 
trum was corrected for atmospheric extinction and flux 
calibrated using observations of spectrophotometric stan- 
dards observed with a l."5 long slit using the same B1200 
grating. A one sigma error array was then extracted ac- 
cording to the FORS reduction pipeline. 

2.2. VIMOS Spectroscopy 

Targets observed with the VIMOS instrument were dis- 
tributed across four slitmasks in the CDF-S, each with a 
field of view of 7' x 8', covering an area of 15' x 16'. These 
observations typically targeted ^5 LAEs and ~40 LBGs 
per mask for a total of 23 LAEs and 164 LBGs. We used 
the VIMOS spectrograph with the MR grism and an OS- 
blue filter mounted at the Nasmyth focus B of UT3 to 
obtain these spectra. The pixel scale of the EEV 4k x 
2k CCD was 0."2 pixel -1 , which, in conjunction with the 
MR grism, led to a dispersion of 2.57A per pixel (R~ 400) 
and a spectral coverage of 4300 - 6800 A on each slitlet. 
The nominal spectral resolution in combination with the 
l."0 slits is 12.5A (FWHM). For these observations, the 
seeing was required to be better than 1". Each mask was 
observed for three hours and eighteen minutes. The VI- 
MOS data were reduced using the standard ESO VIMOS 
pipeline version 2.1.6. For m ore information on t he ob- 
servations and reduction, see ITreister et al.l ([2009). One 
dimensional spectra were then extracted using IRAF's 
apall task. 

The VIMOS spectra were not flux calibrated to a spec- 
trophotometric standard. As a result, we empirically flux 
calibrated the individual spectra based upon the FORS 
spectra. This procedure consisted of fitting a powerlaw 
to the A > 1220A rest-frame continuum of FORS and 
VIMOS spectra via a least-squares algorithm. The in- 
dividual VIMOS galaxy spectral slopes were then ad- 
justed by the difference between the median FORS and 
VIMOS spectral slopes ensuring the same relative spec- 
tral slope. The VIMOS spectra were then normalized 
individually by convolving their fluxes (F>(A)) with the 
V-band transmission curve, then rescaled to match the 
observed V-band magnitude. 

Due to the artificial flux calibration, we also extracted 
the one sigma error spectrum empirically. For each spec- 
trum, a 515A wide top hat, corresponding to 200 pixels 
or 40 resolution elements, was shifted down the spectrum 
wherein, for each pixel, we recorded the median number 
of counts and relative root mean square deviation. Our 
error spectrum is conservative as we only masked the 
Lya emission line and the A5577 skyline, but no other 
features. To check the accuracy of empirically extracting 
the error spectrum, we applied the same routine cover- 
ing the same number of resolution elements to the FORS 
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spectra. This process recovered an error spectrum com- 
parable to the formal error spectrum produced by the 
FORS reduction pipeline. 

2.3. Redshift Identification 

Spectral identification and redshifts were assigned by 
interactively examining the observed spectra in both one 
and two dimensions. All redshifts were independently 
confirmed by two investigators. Because our spectra lack 
systemic redshifts from nebular emission lines, redshifts 
were determined from the Lya emission line. There were 
four continuum-selected galaxies with interstellar absorp- 
tion lines and no Lya emission. For these objects, red- 
shifts were determined from the observed wavelength of 
the Si II 1260A and C II 1334A absorption lines, as they 
are the strongest lines that are not affected by blending. 
We then calculate their expected Lya emission line red- 
shift by applying a median offset of Lya emission and 
interstellar absorption (< v >= 650 km s _1 ) found to 
be characteristic of z— 3 LBGs by S03. For stacking pur- 
poses, we use each object's Lya emission line redshift 
to shift them into the rest-frame. We note that varying 
kinematic offsets between emission and absorption will 
artificially broaden the absorption lines. 

For a large fraction of our LAEs, redshifts are deter- 
mined solely from the Lya emission line. There are many 
dangers involved in single line redshifts, which makes it 
very important to correctly identify this line as high- 
rcdshift Lya. Since we use an [O II] A3727 and an [O 
III] A5007 filter to select redshift z=2.1 and 3.1 LAEs, it 
is unlikely that local [O II] and [O III] emitting galax- 
ies will survive the color cuts. However, intervening 
high-equivalent width [O II] A3727 emission line galax- 
ies (z = 0.34) can be mistaken as z=3.1 LAEs. The [O 
II] doublet AA 3726, 3729 is just resolved in our FORS 
z = 3.1 sample, which would identify the galaxy as a 
low-redshift [O II] emitting interloper. Unfortunately, 
the VIMOS spectral resolution is insufficient to resolve 
the two components of the d oublet, making i t difficu lt to 
distinguish [O II] from Lya. I Gawiser et al.l (|2006aD and 
iGronwall et ail (|2007t) do not find any interlopers among 
LAEs observed with the Inamori Magellan Areal Camera 
and Spectrograph (IMACS) covering a wavelength range 
4000 — 10000A that would show multiple emission lines. 
Nonetheless, as the transmission fr action of the I GM due 
to H I absorption is -0.7 at z = 3 (|Madaulll995l) . we can 
identify [O II] A3727 emitting galaxies from their lack of 
a Lyman decrement. We do not find any low-redshift 
emission line interlopers. Among the targets observed 
with the FORS spectra, we spectroscopically confirmed 
13 of 36 z = 2.1 LAEs, 24 of 55 BXs, 23 of 25 z = 3.1 
LAEs and 3 of 16 LBGs. Of the z~3 VIMOS targets, we 
spectroscopically confirmed 17 of 23 LAEs and 40 of 132 
LBGs. 

3. INDIVIDUAL GALAXY PROPERTIES 

In this section, we analyze the individual galaxy spec- 
tra included in the composites. For individual spectra, 
the faint nature of high-redshift SFGs limits us to study- 
ing trends based on only a few parameters that are mea- 
surable for each galaxy: Lya equivalent width, redshift, 
rest-frame UV color, and rest-UV continuum luminosity. 
For certain galaxies with Lya emission and interstellar 



absorption, we are also able to measure a relative veloc- 
ity offset. Due to the complexities associated with se- 
lection effects, we do not attempt to study the relation- 
ship between rest-UV continuum luminosity and other 
galaxy properties. As a maximum apparent magnitude 
of R<25.5 is used for broadband-selected objects, the 
rest-UV luminosity threshold for z— 2 BX galaxies will 
be ~0.5 magnitudes brighter than for z~3 LBGs. Ad- 
ditionally, we find that 5 BX galaxies, 3 z = 2.1 LAEs, 
and 4 LBGs, i.e., 22% of our R<25.5, Lya emission line 
galaxies, show mu ltiple-peaked Lya emission. Similarly, 
iKulas et all (|2011h found that ~ 30% of their Lya emis- 
sion line galaxies at 2 < z < 3 showed multi-peaked Lya 
emission. We also note seeing evidence for fine-structure 
Si II* emission in the individual spectra of 7 BX galaxies 
and 3 LBGs. 

In 8 of 17 BX galaxies and 15 of 28 LBGs, the only 
visible feature is Lya emission. In 2 BXs and 2 LBGs, 
we observe multiple absorption lines and no Lya emis- 
sion, and in 3 of these 4 galaxies, Lya appears as broad 
absorption. In these cases, interstellar absorption lines 
are used to measure the redshift. In the remaining 7 
BXs and 11 LBGs, Lya emission and interstellar absorp- 
tion are present allowing a measurement of the relative 
velocity offset. For 77 of 81 galaxies, redshifts were mea- 
sured from the Lya emission line. This fraction differs 
significantly from S03, as nearly one third of their LBG 
spectroscopic sample had Lya only in absorption. This 
difference is likely due to lower S/N in our spectra mak- 
ing absorption line-only redshifts difficult. The top left 
panel of Figure Q] shows the distribution of Lya equiv- 
alent widths for our entire spectroscopic sample (black, 
solid line) and UV-bright LAEs (red, dot-dashed line). 
We discuss the impacts of the different Lya equivalent 
width distribution in §6.21 

Among the narrowband-selected galaxies, in all except 
for one, the only visible spectroscopic feature is Lya 
emission. This is likely due to the intrinsically faint na- 
ture of their continua. In the top right panel of Figure 
[U photometric versus spectroscopic rest-frame equiva- 
lent width is plotted for z— 3 LAEs as black stars and 
z~2 LAEs as red diamonds. We find that for z— 3 LAEs, 
photometric equivalent widths are systematically overes- 
timated as the median photometric equivalent width is 
W^phot^ = 69. 2A while the median spectroscopic equiva- 
lent width is Wspoc^ = 36. 6A. However, for z— 2 LAEs, 
these two equivalent width measurements are consistent 
as Wphot,2 = 27. 3A and Ws p0 c,2 = 31. 5A. Photometric 
equivalent widths are calculated by using broadband fil- 
ters (U+B at z—2, and B+V at z— 3) to estimate the con- 
tinuum flux density and a narrowband filter to estimate 
the amount of excess Lya flux. For detai ls on the photo- 
metric equivalent width cal culations, see IGronwall et aLl 
(|2007t) ; iGuaita et ail (|2010D . For particularly faint ob- 
jects, spectroscopic equivalent widths may be biased low 
as a result of fake continuum in individual spectra. Since 
we weight each object in the composites by signal divided 
by noise squared, these objects will have a minimal effect 
on the composite spectra of which all equivalent widths 
are measured spectroscopically. A more detailed analy- 
sis of photometric and spectroscopic equivalent widths is 
necessary to fully understand the cause of this discrep- 
ancy. 
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In galaxy spectra with both Lya emission and interstel- 
lar absorption, we measure the average relative velocity 
offset, Au om _ a bs- In the bottom left panel of Figure[IJ the 
distribution of Av cm _ a b s as a function of Lya equivalent 
width is plotted where the errors represent the range in 
velocity offsets measured from different transitions. The 
median velocity offset is Av cm _ a b s = 584 km s" 1 , similar 
to S03 who found Av cm _ a b s = 650 km s _1 . In all of our 
galaxies, Lya emission is at a higher redshift than in- 
terstellar absorption. This result is consistent with S03 
who find that nearly all of their 323 LBGs show Lya 
at a higher redshift than interstellar absorption. These 
velocity offsets indicate that most likely the emission is 
coming from the farside component of the outflow while 
interstellar absorption, Lya absorption, and blueshifted 
Lya emission are coming from the nearside component. 
Large-scale outflows are caused by feedback from star 
for mation and supernoy ae, as first reported in z^3 LBGs 
bv lSteidel et al.l (|l996a). In general, stellar systemic red- 
shifts are measured from photospheric features from hot 
stars. In the rest-UV, photospheric features are too weak 
to be seen in individual spectra. Systemic redshifts can 
also be determined from rest-optical nebular emission 
lines, but these were not available. Therefore, in the rest 
of our analysis, we measure the relative velocity offset be- 
tween Lya emission and absorption. While velocity offset 
tends to decrease with increasing Lya equivalent width, 
we do not find this trend to be very significant. Further- 
more, none of the UV-bright SFGs with W Lya > 40A 
have visible interstellar absorption lines. 

The spectral slope of the UV continuum is largely 
dependent on the star-formation history and amount 
of extinction. For UV-bright, star-forming galaxies, 
the unreddened UV continuum remains fairly constant 
for ages 10 Myr to 1 Gyr where the spectral slope 
index, 13 (F x oc A* 3 ), r an ges from -2.6 to -2. 1 
(based on models of SFGs from iLeitherer etlrtl Il999f ). 
Both narrowband- and continuum-selected galaxies at 
2 < z < 3.5 have typical ages within t his range 
dPapovich et al.l 120011: IShaplev et al.l 120011: lErb et al.l 
2006bi: lGawiser et al.ll2007HKornei et al.ll2010D indicating 
that large differences in UV spectral slope reflect vary- 
ing amounts of dust extintion. For a large number of our 
galaxies, the individual spectra are of insufficient S/N to 
accurately measure a UV continuum slope. However, us- 
ing broadband photometry, BVRI for z^3 and UBVRI 
for z~2 SFGs, we calculate the UV spectral slope index 
after correcting for IGM extinction and Lya equivalent 
width for each galaxy. The bottom right panel of Figure 
[T] shows the distribution of UV spectral slopes as a func- 
tion of Lya equivalent width for R<25.5, 2 < z < 2.7 
galaxies (blue triangles), R<25.5, 2.7 < z < 3.5 galax- 
ies (red squares), and R> 25.5 galaxies (cyan crosses) 
with a typical err or measuremen t in th e upper right cor- 
ner. LAEs from iGuaita et al.l (120111 ) (black X's) and 
iNilsson et al.1 (|2009f ) (black crosses) are also plotted. We 
find a median spectral slope of (3 — —1.4, then create 
a red ((3 > —1.4) and a blue (f3 < —1.4) spectral slope 
subsamplc. Galaxies with nondctcctions in the B, V, 
and/or R bands have been omitted from this plot. We 
find a large range in (3 values from f3 = —2.7 to f3 = 0.1, 
however, for galaxies with rest-frame Wl Y q > 20A, j3 
ranges from j3 = —2.7 to j3 = —0.8 indicating a rela- 



tion between stronger Lya emission and bluer spectral 
slopes. The minimum observed spectral slopes are in 
good agreement with model predictions for a dust-free 

SFG with an age of 10 Myr (j3 2.6, ILeitherer et~aTl 

119991 ). UV-bright SFGs with W LyQ < 20A subset have a 
mean spectral slope of < (3 > = —1.22 ± 0.47 while UV- 
bright SFGs WLy Q > 20A have a bluer mean spectral 
slope of < (3 >— —1.52 ±0.45 where errors represent the 
standard deviation on the mean. Moreover, the red and 
blue spectral slope subsamples have median Lya equiva- 
lent widths of W = 15. lA and W = 32. 2A respectively. 
These trends display a relation between increasing Lya 
equivalent width and bluer spectral slope, which has pre- 
viously been observed in z~3 LBGs (IShaplev et al.ll2~003l : 
ISteidel et al1l2010t iKomei et al.ll2010D 

Through dividing the sample into two redshift bins, 
2.0 < z < 2.7 and 2.7 < z < 3.5, we find that galax- 
ies at both redshifts show a similar relationship between 
Lya equivalent width and spectral slope. Surprisingly, 
galaxies at z~2 show bluer continuum slopes on av- 
erage (< (3 >= —1.87 ± 0.51) than galaxies at z^3 
(< (3 >= — 1.34 ± 0.33), which appears as a systematic 
offset independent of Lya equivalent width and is further 
discussed in § 16.51 This trend is still observed when us- 
ing the same observed photometry (BVRI) to determine 
rest-UV slopes at both redshifts. 

4. STACKING PROCEDURE 

In generating the following composite spectra, we be- 
gin by identifying a subsample of galaxy spectra to be 
combined. We exclude LAEs whose spectra are plagued 
by poor sky subtraction or have minimal wavelength cov- 
erage. Of the 52 spectroscopically confirmed LAEs, we 
include 26 observed with the FORS instrument and 10 
with the VIMOS instrument. Additionally, SFGs lack- 
ing precise redshifts determined by either the Lya emis- 
sion line or several interstellar absorption lines are also 
excluded from the composite. These criteria yield 28 
VIMOS LBGs and 17 FORS BXs out of 73 spectroscopi- 
cally confirmed continuum-selected galaxies. The effects 
of omitting LBGs without detectable Lya emission red- 
shifts is discussed in section [6. II 

We proceed to shift the extracted, one-dimensional, 
flux-calibrated spectra into the rest frame. The spectra 
are subsequently rebinncd to a dispersion of lA per pixel 
corresponding to roughly the FORS resolution. We then 
renormalize them to their 1250 — 1300A continuum. Out- 
lying data points due to cosmic-ray events and sky sub- 
traction residuals are systematically masked, totaling less 
than 0.2% of the data. Finally, using an optimal weight- 
ing routine where weights are defined as 1250 - 1300A 
continuum strength divided by noise squared, w e co- 
add the individual spectra ()Gawiser et al.l 1200 6a) . We 
note that a variety of UV spectral slopes in conjunction 
with the majority of the VIMOS spectra covering only 
A < 1600 A causes the S/N to drop rapidly and the con- 
tinuum level to become uncertain above this wavelength. 

In order to understand the connection between UV- 
bright star-forming galaxies with and without strong 
Lya emission, we divided our dataset into three sub- 
samples: UV-bright non-LAE, UV-bright LAE, and UV- 
faint LAE. In Table [1] we present our selection cri- 
teria for and the number of galaxies in each com- 
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posite. We also define a UV-bright SFG subsample 
with R<25.5 which is composed of the UV-bright non- 
LAE and UV-bright LAE subsamples. The UV-bright 
non-LAE and UV-bright LAE subsets differ slightly in 
magnitude, < Ruv-brightnon-LAE >= 24.46 and < 
Ruv-brightLAE >= 24.84 with a scatter of ±0.68 and 
±0.42 respectively, although this difference is not signifi- 
cant. This difference is in part due to the fact that larger 
Lya line fluxes make spectroscopic identification easier. 
Selection effects based on luminosity and Lya equivalent 
widths are discussed in detail in §6.11 The primary goal 
for these different spectroscopic subsamples is to under- 
stand the effects of increasing Lya emission and decreas- 
ing magnitude on galaxy properties. We do note that 
our R< 25.5 LAE subsample contains galaxies with a 
broad range of properties including galaxies with strong 
interstellar absorption lines, Si II* fine-structure emis- 
sion, and faint continua where the only feature seen is 
Lya emission. We also subdivide our dataset into a red 
(/3 > —1.4) and blue (j3 < —1.4) spectral slope subsam- 
ple based on their intrinsic rest-UV colors. We exclude 6 
narrowband-selected galaxies with B, V, or R band mag- 
nitudes greater than 1 magnitude from these composites 
as their spectral slopes are not well constrained. 

5. ANALYSIS OF COMPOSITE SPECTRA 

The absorption and emission features observed in these 
spectra are produced in gas associated with galactic- 
scale outflows, H II regions, and the winds and photo- 
spheres of massive stars. Figure [5] shows the UV-bright 
SFG composite spectrum of 59 R<25.5 SFGs. The most 
prominent features in the rest-UV composite spectra are 
blueshifted interstellar absorption and H I Lya. These 
features trace large-scale outflows of interstellar medium, 
which have been modelled in detail for high-redshift 
LAE s and LBGs (|Verhamme et all 120081 : lLaursen et al.l 
1201 If) . In individual galaxy spectra, due to low S/N, only 
redshifts can be determined from strong interstellar ab- 
sorption. However, with higher S/N, composite spectra 
allow for more robust measurements of their equivalent 
widths as well as revealing fainter features such as neb- 
ular emission. 

In Figure [31 the top panel displays our composite UV- 
bright non-LAE spectrum (27 UV-bright non-LAEs), the 
middle panel shows our UV-bright LAE composite spec- 
trum (32 UV-bright LAEs) and the bottom panel shows 
our UV- faint LAE composite spectrum (22 LAEs). For 
our UV-faint LAE composite spectrum, limited spectral 
coverage in conjunction with all of these LAEs being at 
z = 3.1 causes our sample size to drop to 8 galaxies 
(35%) at A < 1340A. In addition, we have masked the 
region around 1360A due to residuals from the A5577A 
skyline. For this reason, we limit our analysis of the LAE 
composite to wavelengths blueward of 1340A. 

Across the UV-bright LAE and UV-bright non-LAE 
composites, the most striking trends are the decrease 
in interstellar absorption and increase in Si II* nebu- 
lar emission with increasing Lya emission. In this sec- 
tion, we first focus our analysis on the Lya emission and 
absorption line profile. We then discuss low- and high- 
ionization interstellar absorption line features. Finally, 
we proceed to discuss nebular features and fine-structure 
Si II* emission within the systemic component of the 



galaxies. 

5.1. Lyman Alpha Emission and Absorption 

Lya photons originate from recombinations in H II 
regions around massive stars, yet a number of effects 
impact the strength, shape, and relative velocity of the 
emergent Lya profile; including intrinsic star formation 
rate, amount of dust, dust geometry, covering fraction, 
and the kinematics of the outflowing gas. For a more 
in depth discussion on the shape of the Lya profiles in 
this dataset, see Gawiser et al. (2012) [in preparation]. 
Nonetheless, we are still equiped with the overall strength 
of Lya emission and absorption as well as its relative ve- 
locity offset with respect to other spectroscopic features. 

For our UV-bright SFG composite (Fig. [2]), we mea- 
sure a rest-frame Lya emission equivalent width of W Q = 
19.1 ± 1.5A. We also note that there is evidence for weak 
blueshifted emission at a velocity offset of Aw cm _ a b s ~ 
5000 km s -1 separated by Lya absorption. Lya photons 
are much more likely to escape from a galaxy if their 
relative velocities have been shifted off resonance with 
respect to the bulk component of the neutral material. 
We find it interesting that there is blueshifted emission 
at a very high velocity offset as the UV-bright SFG com- 
posite represents an average of all the UV-bright galax- 
ies including several with strong Lya absorption. For the 
main Lya absorption, we measure an equivalent width of 
W = -2.7±0.3A at a velocity offset of Aw cm _ abs ~ 2400 
km s _1 . Equivalent widths and velocity offsets for our 
UV-bright SFG composite spectra are presented in Table 
HJ Since the composite spectra include galaxies observed 
with the FORS and VIMOS instruments, we smooth all 
spectra by the VIMOS resolution (v ~ 650 km s _1 ), 
which we calculate to be 3.0 A under optimal observ- 
ing conditions. We also set an upper bound of 3.1 A 
on our spectral resolution as the minimum FWHM value 
measured among the strong interstellar features. Assum- 
ing an effective spectral resolution of 3.0 A and subtract- 
ing the instrument FWHM from the observed FWHM 
in quadrature, we compute a deconvolved FWHM(Lya) 
= 530 ± 30 km s^ 1 in the UV-bright SFG composite. 
While the UV-bright SFG composite represents the av- 
erage strength of Lya emission among our R<25.5 SFGs, 
this is not representative of the entire population of SFGs 
due to selection effects. For instance, spectroscopic iden- 
tification is easier using an emission line rather than an 
absorption line. As objects become fainter, the fraction 
of galaxies with Lya in emission increases. We discuss 
the implications of this effect in §53] Additionally, the 
weight causes brighter galaxies to contribute more to the 
composite. These factors have opposing effects on the 
strength of the Lya feature in Figure [21 so it will not 
necessarily be that of the average SFG. However, we do 
find the median Lya equivalent width of the UV-bright 
SFG sample (Wo = 2 1.0 A) to be consistent with that 
of the UV-bright SFG composite. A wide distribution 
of Lya equivalent widths and line profiles is seen in our 
individual spectra, which we discuss in section 16.21 In 
§6.3[ we discuss the relationships among the strength of 
the Lya line and other galaxy properties. 

In the UV-bright non-LAE composite (top panel of 
Fig- El): we observe Lya both in emission and blueshifted 
absorption with an average velocity offset of Au cm _ a b s ~ 
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2400 km s 1 . The blueshifted absorption extends from 
the Lya emission line to Ai> cm - a bs > 4500 km s -1 with 
an equivalent width of W D = —4.9 ± 0.9A. For the emis- 
sion component, we measure a Lya emission line equiv- 
alent width of W = 8.4 ± 0.9A, which we find to be 
comparable to the second highest Lya equivalent width 
quartile of S03 (W = 11.00 ± 0.7lA). Our results for 
the UV-bright non-LAE and UV-bright LAE composites 
are presented in Table [3] The observed FWHM is right 
at our instrument resolution, so we set a 3<r upper limit 
of FWHM(Lya) < 200 km s" 1 . 

The UV-bright LAE composite (middle panel of Fig. [3]) 
has a complex Lya profile, as it includes galaxies with 
rcdshifted absorption and multiple-peaked emission. For 
the main Lya emission line component, we measure de- 
convolved FWHM (Lya) = 410 ± 100 km s" 1 and an 
equivalent width of W D = 50.0 ± 4.8A. The faint 
blueshifted emission at very high velocity offset in this 
composite is the origin of the same feature in the UV- 
bright SFG composite. We measure a velocity offset of 
Aw cm _ a bs ~ 4600 km s , which agrees with the maxi- 
mum velocity offset of the Lya absorption in the UV- 
bright non-LAE subsample. In the individual galaxy 
spectra contributing to the UV-bright LAE composite, 
we do not find any significant emission at this velocity off- 
set, but we do note that the stacking procedure can reveal 
features too faint to be seen in individual spectra. The 
majority of Lya photons will be absorbed at the relative 
velocity of the bulk of the outflowing neutral medium. 
However, Lya photons farther in velocity s pace from the 
line ce nter are less likely to be re-absorbed. iSteidel et al.1 
(2010) find outflows around z = 2.3 SFGs to best be char- 
acterized by a powerlaw with larger velocities at larger 
radii. If the blueshifted emission at Av om _ a b s ~ 4600 km 
s _1 is real, it reveals extended e mission of Lya phot ons 
from a diffuse component of H I. ISteidel et al.l (|2011l ) re- 
cently found extended Lya emission in a sample of 92 
continuum-selected galaxies at z~2.65. 

In our UV-faint LAE composite (bottom panel of 
Fig- El, we see no strong evidence for blueshifted emis- 
sion. We measure a Lya equivalent width of Wq = 
42.1 ± 6. 2 A and a velocity FWHM of FWHM (Lya) = 
380 ± 110 km s -1 , consistent with the UV-bright LAE 
composite. 

5.2. Low-Ionization Absorption Lines 

Blueshifted singly-ionized metal absorption lines probe 
the neutral component of outflowing interstellar mate- 
rial on the nearside of the galaxy. Among our compos- 
ite spectra, there are a host of low-ionization interstellar 
absorption lines. Due to limited S/N, we restrict our 
analysis to the strongest low-ionization lines that we de- 
tect at high significance: Si II A1260, O I A1302 + Si II 
A1304, C II A1334, Si II A1526. All of these features have 
been well-s t udied in high-redshif t LBGs dPettini et al.l 



2008 



2010; 



2000 . 120021: iShaplev et al.l 120031: ICabanac et al.l 
Quider et al.l l2009t iDessauees-Zavadskv et al.l _ 
Talia et al.ll2011l:lJones et al.ll201 ]) and nearby starburs t 
galaxies (Hcckma n et al.l Il998t iLeitherer et al.l 1201 lh . 
We also find Fe II A1608 and Al II A1670 at low sig- 
nificance, but they are located in a region of the spec- 
trum where we have fewer galaxies contributing to the 
composite and therefore lower S/N. For these reasons, 



we only record their equivalent widths. Table [3] reports 
measurements of equivalent widths and average velocity 
offsets of the interstellar lines. As we observe a distri- 
bution of velocity offsets among our individual galaxies, 
the FWHMs of all transitions other than Lya will be 
artificially broadened to some extent. Due to the intrin- 
sically faint nature of our UV-faint LAE galaxies' con- 
tinua in combination with continuum normalization, we 
do not observe any significant interstellar features in the 
UV-faint LAE composite. As a result, we restrict our 
analysis to the UV-bright LAE and UV-bright non-LAE 
composites for the rest of the analysis. 

The strong low-ionization interstellar lines are useful 
for measuring the kinematics of the neutral gas compo- 
nent within the outflow. Figure [4] and Figure [5] show 
zoomed in regions around the strong absorption lines 
for the UV-bright LAE composite (top panel) and UV- 
bright non-LAE composite (bottom panel). The ver- 
tical lines indicate the rest-frame wavelength of each 
transition with respect to Lya emission. Relative to 
Lya, we measure an average velocity offset of the low- 
ionization features of Av cm - a bs = 640 ± 10 km s" 1 and 
Aw om _ ab s = 610 ±60 km s" 1 for the UV-bright non-LAE 
and UV-bright LAE composites. Using the Si II A1260 
and C II A1334 absorption lines, which are the least af- 
fected by blending, and assuming a spectral resolution of 
3.0A, we measure a deconvolved FWHM = 430±60 km 
s _1 for the UV-bright SFG composite. We subsequently 
measure a FWHM in the UV-bright LAE and UV-bright 
non-LAE composites of FWHM = 420 ± 80 km s^ 1 and 
FWHM = 340 ± 100 km s" 1 respectively. The differ- 
ent velocity offsets in the UV-bright LAE and non-LAE 
composites is likely biasing the FWHM in the UV-bright 
SFG composite to larger values. Our results are con- 
sistent with S03, who stacked 811 LBGs at a spectral 
resolution of 2.6A and found an average FWHM(LIS) = 
450 ± 150 km s" 1 . 

We measure the degree of saturation of the Si II 
transitions by comparing the equivalent widths of Si 
II A1260 and Si II A1526. On the linear part of the 
curve of growth, the ratio of W o (1260)/W o (1526) > 
5, while we measure W (1260)/W o (1526) = 1.1 ± 
0.2 in the UV-bright SFG composite, consistent with 
unity and the material being optically thick. This ra- 
tio is W o (1260)/W o (1526) = 2.4 ± 1.5 in the UV- 
bright LAE composite, as the A1526 transition is 
marg inally dete c ted, ma king t he uncertainty quite large. 
S03, lErb et alJ (|2006bD . and ISteidel et al. l (120101) find 
the Si II transitions to be consistent with satura- 
tion in each of their composites, and S03 report a 
W o (1260)/W o (1526) = 0.95 in their all-LBG composite. 
Since the strong low-ionization transitions are all satu- 
rated, they are not useful for determining chemical abun- 
dances. Metal abundances can be derived from weaker 
features lying on the linear part of the curv e of growth in 
conduction with H I column densities (e.g., iPettini et al.1 
2000, 2002). Wc do not observe any of these features with 
high enough S/N and cannot measure H I column den- 
sities, so we foc us our analysis on the strong transitions. 
We do note that lErb et al.( (|2006af) used the metallic ab- 
sorption lines near A1370 and A1425 to confirm their rest- 
frame optical metallicity measurements. We find these 
regions in our UV-bright SFG composite to be consis- 
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tent with the low- mass composite of lErb et al.l (12 0068) 
which had a metallicity of Z < 0.33Z ? . iFinkelstein et al.1 
( 2009) measured limits on the metallicities of two LAEs 
using rest- frame optical emission lines to be Z < 0.1 7Zq 
and Z < 0.28Z Q . 

5.3. High- Ionization Absorption Lines 

In addition to low-ionization absorption, we also de- 
tect at high significance several high-ionization interstel- 
lar absorption lines, namely N V AA1238, 1242, Si IV 
AA1393, 1402, and C IV AA1548, 1550. These transitions 
trace gas at T > 10 K that has been ionized by feedback 
from supernovae, stellar winds, and collisional processes 
associated with the outflow. They are therefore useful 
for probing the ionized component and ionization state 
of the outflow. Among all of our composites, we find that 
the high-ionization line profiles are dominated by the in- 
terstellar component with a smaller contribution from 
the stellar component. We analyze the most prominent 
high-ionization features first (Si IV and C IV), then in- 
terpret the N V line profile, which we detect at lower 
significance. 

In the UV-bright non-LAE composite, for the AA 
1393, 1402 Si IV transitions, we measure comparable ve- 
locity offsets of Au om _ a bs ~ 560 km s" 1 and Au cm _ a b s ~ 
590 km s _1 . In contrast, for the UV-bright LAE com- 
posite, we measure velocity offsets of Ai; om _ a b s = 440 
km s" 1 and Au em _ a b s = 870 km s _1 respectively. The 
UV-bright LAE composite shows strong absorption at 
A1393 while having marginally detected absorption at 
A1402 likely due to redshifted emission from a P Cygni 
component filling it in, and thereby biasing our veloc- 
ity offset to a higher value. We also find the signature 
of a P Cygni profile in the UV-bright non-LAE compos- 
ite, although at lower significance and it does not ap- 
pear to have a significant effect on the velocity offset of 
the A1402 absorption. For a doublet ratio of 2:1, the Si 
IV AA 1393, 1402 transition is on the linear part of the 
curve of growth and therefore optically thin. In the UV- 
bright SFG composite, this ratio is 1.9:1, however, this 
ratio is 3.5:1 in the UV-bright LAE composite and 1.4:1 
in the UV-bright non-LAE composite. The higher dou- 
blet ratio in the UV-bright LAE composite is again due 
to P Cygni emission filling in A1402 absorption. In the 
UV-bright non-LAE composite, continuum uncertainty 
due to broad stellar absorption makes measuring this ra- 
tio difficult and may bias this ratio towards 1:1. Thus, 
we conclude that the Si IV transition appears consistent 
with being optically thin in all composites. 

For the UV-bright SFG composite, we find the rela- 
tive velocity offset of the Si IV AA 1393, 1402 transi- 
tions (A^em-abs ~ 590 km s _1 ) to be in good agree- 
ment with the average velocity shift of the low-ionization 
absorption lines (A?; em _ a bs = 550 km s _1 ). Moreover, 
the average deconvolved FWHM of the Si IV doublet is 
FWHMsi/y = 420 i 110 km s" 1 , also consistent with 
that of the low-ionization absorption lines FWHMi/s = 
360 ± 60 km s" 1 . 

In the UV-bright SFG composite, for C IV blended 
at A1549.5, assuming a doublet ratio of 1:1, we mea- 
sure Av cm _ a bs = 800 km s^ 1 , significantly larger than 
the low-ionization lines. However, the velocity width of 
FWrlMc/y = 450 ±60 km s _1 is consistent with the low- 
ionization lines. There appears to be some broadening 



due to stellar winds along with redshifted emission char- 
acteristic of a P Cygni profile, which may be biasing the 
velocity offset to larger values. For both the UV-bright 
LAE and UV-bright non-LAE composites, we find com- 
parable velocity offsets to the UV-bright SFG composite, 
which are in all cases larger than that of the Si IV A1393 
transition and the low-ionization absorption lines. 

In the UV-bright LAE composite, we see weak C IV 
absorption, yet in the UV-bright non-LAE composite, we 
find significantly stronger C IV absorption corresponding 
to an increase in absorption strength of ~100%. S03 also 
report seeing more C IV absorption in their WL yQ = 11 A 
composite than their Wl Y q = 52A composite, although 
they only measure an increase in strength of ^30%. The 
relatively weak C IV absorption seen in the UV-bright 
LAE composite is particularly interesting given that the 
strength of the Si IV A1393 absorption is comparable in 
both the UV-bright LAE and UV-bright non-LAE com- 
posites, and the velocity FWHMs of the high-ionization 
transitions are comparable across all composites. We fur- 
ther note from the relative strength of the interstellar 
C II A1334 absorption to C IV A1549 absorption that 
the outflowing component of carbon appears to have a 
higher covering fraction in a lower ionization state, as it 
is unlikely that the velocity FWHM of the neutral gas is 
higher than that of the ionized gas. 

Significant N V absorption is visible in the UV-bright 
SFG composite with a large velocity FWHM (FWHM~ 
1200 km s _1 ) and velocity offset (Av em _ a b s = 1150 km 
s _1 ). The N V line profile shape is consistent with that 
of a P Cygni profile, indicating that there is both a stellar 
and an interstellar component. The velocity offset of the 
N V absorption is significantly larger in the UV-bright 
non-LAE composite (Av om _ a b s = 1300 km s" 1 ) than the 
UV-bright LAE composite (Av cm _ a b s = 600 km s _1 ). 
In spite of this difference in velocity offset, the velocity 
width is consistent among all three composites (FWHM~ 
1200 km s _1 ). We interpret the larger velocity offset 
in the UV-bright non-LAE composite as evidence for a 
significant amount of interstellar N V absorption in the 
outflowing component. In this composite, both the N 
V velocity offset and velocity FWHM are significantly 
larger than that of the Si IV and C IV absorption. 

5.4. Stellar Features 

The radiation pressure from hot m assive stars can gen- 
erate winds of 2000-3000 km s _1 (IGroenewegen et aLI 
1989). These winds are seen spectroscopically as broad 
absorption in low density environments or P Cygni pro- 
files, blueshifted absorption followed by redshifted emis- 
sion, in higher density environments (|Leitherer et al.l 
1995). In the far UV, the most prominent stellar fea- 
tures are N V AA1238, 1242, Si IV AA 1393, 1402, C IV 
AA 1548, 1550, and He II A1640. There is no evidence for 
He II A1640 emission in any of the composites, however, 
the S/N in that region is much lower due to limited spec- 
tral coverage and an uncertain continuum level. For the 
Si IV and C IV transitions, a combination of strong in- 
terstellar and photospheric absorption contribute to the 
line profiles, which can be difficult to differentiate. In 
the UV-bright LAE composite, we find evidence for O 
III] A1664 nebular emission, which is not detected in the 
UV-bright non-LAE composite. Due to low S/N and an 



Spectroscopy of High-Redshift Galaxies 



9 



uncertain continuum level, we do not attempt any further 
analysis. This result is qualitatively consistent with S03 
who reported significantly stronger O III] A 1664 emission 
in their highest Lya equivalent width quartile relative 
to their lower three Lya equivalent width quartiles. A 
number of weak stellar features are also visible in the UV- 
bright SFG composite, the strongest being C III A1176 
absorption at a velocity offset of Av cm _ a b s = 380 km s _1 , 
which further confirms that Lya is typically redshifted 
from the systemic redshift of SFGs at thes e redshifts by 
~400 km s _1 in excellent agreement with ISteidel et all 

The stellar wind component of the Si IV transition is 
only significant for blue giant and supergiant stars. In the 
UV-bright SFG composite, we observe extended absorp- 
tion that ranges from zero velocity offset to Au cm _ a b s = 
1600 km s _1 as well as redshifted emission characteris- 
tic of a P Cygni profile, consistent with high velocity 
stellar winds from massive stars. We observe a simi- 
lar profile in the UV-bright LAE composite except with 
redshifted emission filling in the Si IV A1402 absorption. 
Although the UV-bright non-LAE composite shows weak 
redshifted emission, there does not appear to be signif- 
icant broad absorption. Given the strength of Si IV A 
1393, the most significant difference between these two 
composites is the much lower amount of Si IV A 1402 in- 
terstellar absorption in the UV-bright LAE composite. 
The shape and strength of these line profiles indicates 
the presence of a population of young blue giant and su- 
pergiant stars. 

The stellar wind contribution to the C IV profile is 
most dependent on the presence of main-sequence, giant 
and supergiant O stars. We observe strong C IV AA 1548, 
1550 absorption in the UV-bright non-LAE composite, 
though it is difficult to disentangle the stellar wind com- 
ponent from the interstellar absorption component. We 
find evidence for weak extended absorption in the UV- 
bright SFG and UV-bright non-LAE composites, but not 
in the UV-bright LAE composite. We also find significant 
redshifted emission in the UV-bright non-LAE compos- 
ite. A P Cygni profile i ndica t es the presence of M >30M^ 
stars (jLeitherer et all fT995l iPettini et all 120001 ). This 
emission is also seen in the UV-bright SFG composite. 

In both the UV-bright LAE and UV-bright SFG com- 
posites, there is significant blucshiftcd absorption and 
redshifted emission in the N V line profile, characteristic 
of a P Cygni profile. We interpret the low velocity off- 
set (Aw cm _ a b s = 600 km s _1 ) and large velocity FWHM 
(FWHM~ 1200 km s" 1 ) of N V absorption in the UV- 
bright LAE composite as evidence for stellar emission 
an d absorption dominating the N V profile (see Fig. 1 
in iGroenewegen et al.| [l989). Given the shape of the N 
V line profile in the UV-bright SFG composite, we in- 
terpret this as evidence for a significant amount of high- 
ionization absorption in the outflow as well as emission 
and absorption due to stellar features at the systemic 
redshift. 

5.5. Si II* Emission Lines 

An advantage of composite spectra is the ability to 
reveal weak features that are typically undetected in in- 
dividual spectra. Among these features, the most promi- 
nent are the Si II* fine structure emission lines at A 1265, 
A1309, and A1533. These transitions have also been ob- 



served in z~3 LBGs by S03, z~4 LBGs by IJoiies et al.l 
(|20T1 . and in a local LBG bv lFrance etldl (|2010l ). 

The close proximity of these lines to strong absorp- 
tion features (Si II A1260, O I + Si II A1303, and Si II 
A1526) may attenuate the blue edges, biasing the cen- 
troids to the red. Our results for Si II* emission in the 
UV-bright SFG composite are presented in Table |U In 
the UV-bright non-LAE composite, we observe signifi- 
cant Si II* A 1265 emission yet we do not see the A 1309 
and A1533 transitions. For the A1265 emission, we mea- 
sure an equivalent width of Wo = 0.5 ± 0.2A at a veloc- 
ity offset of Au cm _ a bs = 300 km s _1 . We find slightly 
stronger Si II* A1265 in the UV-bright LAE composite 
with an equivalent width of Wo = 0.7 ± 0.2 A at a ve- 
locity offset of A?; em _ a b s = 120 km s . In this compos- 
ite, the Si II* A1265 emission is double-peaked, which is 
likely due to only individual galaxies having Si II* emis- 
sion at different velocity offsets. We observe A1533 emis- 
sion at a much higher signficancc in the UV-bright LAE 
composite than the UV-bright non-LAE composite with 
an equivalent width of Wo = 1.2 ± 0.4 A compared to 
Wo < 0.3A at velocity offsets of Aw cm _ a b s = SOkms" 1 
and Aw cm _ a bs = 380kr7is _1 respectively. We do not find 
any evidence for A 1309 emission in a ny of the composite 
spectra. In contrast, both S03 and Uones et al.l ([20111 ) 
observe Si II* A 1309 emission in their composite spectra. 

Si II* emission likely stems from photoionized regions 
around massive stars in the systemic component of the 
galaxy where electron densities are N e = 10 2 — 10 3 cm~ 3 
with temperatures of T = 10 4 K. In this environment, the 
Si II and Si III relative abundances are determined from 
radiative recombinations for which the recombination 
rates are comparable to the collisional excitation rates 
(jShull fc van Steenberdll982fl . S03 find, through model- 
ing the LBG nebular emission lines with the CLOUDY96 
software package, that all models that simultaneously 
match their O, C, and H/3 line ratios overpredict the Si 
II emission-line strengths by more than an order of mag- 
nitude. This result leads them to conclude that the Si 
II emission lines cannot originate from H II regions. As 
we do not observe these other transitions, we are unable 
to comment on the plausibility of Si II* fine structure 
emission stemming from H II regions. 

A different explanation for Si II* emission is that these 
lines originate instead from the outflowing component. 
In this scenario, each UV photon absorbed at a resonance 
transition (A = 1260, 1304, 1526) or a fine-structure tran- 
sition must be re-emitted through either a resonance 
or a fine-structure transition. In the absence of dust, 
the sum of the resonance and fine-structure emission 
should equal the amount of absorption. Our UV-bright 
non-LAE composite has significantly stronger resonance 
absorption than fine-structure emission while our UV- 
bright LAE composite shows a much larger ratio of fine- 
structure emission to resonant absorption. In contrast, 
S03 observe the fine-structure emission lines to be an or- 
der of magnitude weaker than the resonance lines. The 
discrepancy between the emission and absorption compo- 
nents could be due to the presence of dust extinguishing 
Si II photons during resonant scattering. Furthermore, 
we observe the Si II* fine-structure lines to have narrow 
velocity widths, whereas we would expect the emission 
to come from the full range of outflow velocities and to 
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therefore be as broad as the Lya emission line. For these 
reasons, we find it unlikely that the Si II fine -structure 
lines o riginate from the outflowing component. lErb et all 
(2010) also discuss possible mechanisms for the origin of 
fine-structure emission lines, but find no satisfactory ex- 
planation. 

6. DISCUSSION 

Combining the spectra of both continuum- and 
narrowband-selected galaxies has allowed us to study the 
average properties of 2 < z < 3.5 SFGs. Among these 
properties, we are able to examine trends based on Lya 
equivalent width, rest-frame UV color corrected for IGM 
extinction, redshift, and R band luminosity for each in- 
dividual galaxy. Due to the complexities associated with 
selection effects, we do not attempt to draw any infer- 
ences from different R band luminosities. We have di- 
vided our sample based upon Lya equivalent width for 
galaxies with R<25.5 and will first examine trends as- 
sociated with different amounts of Lya emission. Due 
to the faint nature of our objects' continua, we parame- 
terize their observed colors in terms of their photometri- 
cally measured rest-UV spectral slope after correcting for 
IGM extinction and Lya equivalent width, then divide 
our sample in half based on this parameter. Additionally, 
our spectroscopic sample divides evenly into two redshift 
bins (2 < z < 2.7, 2.7 < z < 3.5) allowing us to look for 
any trends across this redshift interval. This allows us 
to study trends in spectroscopic properties in relation to 
spectral slope and redshift. 

6.1. Selection Effects 

As our subsamples of galaxies are determined on the 
basis of different selection criteria, it is important to iden- 
tify any photometric or spectroscopic biases that may 
affect our analysis of the underlying galaxy populations. 
For continuum-selected galaxies, one of the most influ- 
ential effects is the dependence of color selection criteria 
on luminosity and redshift due to Lya forest absorption. 
Another major bias stems from the prevalence of spec- 
troscopic redshifts determined solely from Lya emission 
in conjunction to a minimum Lya flux necessary for Lya 
emission to be seen. For galaxies without distinct ab- 
sorption lines, this causes spectroscopic identification to 
be Lya flux dependent. 

The selection process for our continuum-selected ob- 
jects creates two major biases that affect our resultant 
galaxy population. First of all, our continuum-selected 
LBGs were chosen to have V-R < 1.2 thereby target- 
ing objects within a redshift range of 2.6 < z < 3.5. 
However, the V band is extinguished by the Lya forest 
making it redshift dependent. At z — 2.5, the amount 
of IGM extinction in the V band is AV= 0.0 mag, while 
at z — 3.5, AV= 0.5 mag. For this reason, the range of 
intrinsic colors of the continuum-selected population is a 
function of redshift. Secondly, the colors of continuum- 
selected galaxies is also a function of R band magnitude. 
The combination of requiring our targets to have very 
faint U band magnitudes (([/ — V)ab > 1-2) and U band 
images having a finite depth will indirectly restrict the 
range of V-R color to bluer colors at fainter luminosities. 

Spectroscopic success rate is not only dependent on ap- 
parent magnitude, but also on the fact that it is easier 
to measure a redshift from an emission line rather than 



an absorption line. Our spectroscopic incompleteness as 
a function of R band magnitude is plotted in Fig [6] for 
narrowband-selected objects (top panel) and broadband- 
selected objects (bottom panel). In each panel, the black 
line indicates the total number of objects observed spec- 
troscopically and the red dot dashed line shows the num- 
ber of spectroscopically confirmed galaxies. The blue 
dashed line shows the number of spectroscopically con- 
firmed galaxies without Lya emission. Our confirmation 
rates are 63% for narrowband-selected objects and 33% 
for broadband-selected objects. Only 5% of broadband- 
selected objects have redshifts determined solely from 
absorption lines. Due to the low S/N of our individual 
spectra, we found it difficult to robustly measure red- 
shifts for galaxies without Lya emission. Galaxies with- 
out Lya emission line redshifts likely have Lya equivalent 
widths below a certain threshold as a function of lumi- 
nosity. We assume thi s to b e consistent with no Lya 
emission. ISteidel et all pOOOf ) find that roughly half of 
their z^3 spectroscopically-confirmed LBGs show Lya in 
absorption. Assuming a similar distribution among our 
broadband-selected objects, this would imply a spectro- 
scopic confirmation rate of 56%. 

6.2. Individual Galaxy Lya Profiles 

Recently, lLaursen et al.l (|201lD used cosmological hy- 
drosimulations to examine the effect of the IGM on 
the Lya line profile escaping high- redshift, star- forming 
galaxies. Their findings qualitatively agree on the shape 
and strength of the redshifted Lya emission line, how- 
ever, they predict more blueshifted emission than has 
typically been observed for high-redshift SFGs. At a 
redshift of z — 3.5, they find the effect of the IGM 
is insufficient to explain the lack of a blueshifted Lya 
emission line. For our sample of 12 multiple-peaked Lya 
emission line galaxies, we find a mean velocity offset of 
< Av ve d-uue >= 660 ± 2 29 km s -1 . This m easurement 
is in good agreement with lKulas et al.l (|201lD who found 
a mean velocity offset of < Av ie d-biue > = 741 ± 39 
km s _1 (errors represent the standard deviation of the 
mean). Ho wever, we do find a larger range of velocity 
offsets than iKulas et al.l (|201lD from Av le d-biue = 330 
km s _1 to Av lec i-biue — 1080 km s" 1 . We note that 
this result is consistent with ISteidel et al.l (|2010| ) who 
found a broad range of Lya emission line velocity offsets 
(Av = +485 ±175 km s _1 ) with respect to the systemic 
redshift. This distribution is also comparable to the ob- 
served velocity offset between Lya emission and inter- 
stellar absorption seen in individual spectra (Fig. [TJ . 

Among the spectroscopic sample, there is a broad dis- 
tribution of Lya line strengths and profile types, which 
fall into four categories: emission only, absorption only, 
both emission and absorption, and neither emission nor 
absorption. Out of the continuum-selected galaxies in 
our composites, we find 45% have Lya emission lines 
strong enough to be classified as LAEs (WL ya > 20A) 
while only 10% show Lya only in absorption. This 
fraction of UV-bright LAEs is roughly a factor of two 
larger than the 20 — 25% that S03 report for their z^3 
LBGs. However, due to the poor S/N of our data, it 
was difficult to determine redshifts for LBGs without 
Lya in emission and we remove any galaxies lacking 
precise redshifts. Removing all galaxies without Lya 
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emission, our fraction of LBGs that would be classified 
as L AEs rises to 5 0% a nd the fraction of LBG-LAEs 
from iSteidel et all (|2000D increases to -40%. Exclud- 
ing LBGs with Lya in absorption, we report a median 
Lya equivalent width of W^ya = 20. OA. Similarly, two 
of our R<25.5 narrowband-selected objects did not yield 
spectroscopic Lya equivalent widths high enough to be 
classified as LAEs and were included in the UV-bright 
non-LAE stack. We find that 24% of our continuum- 
selected galaxies have both Lya in emission and absorp- 
tion. For these objects, we use their emission line equiv- 
alent widths in our analysis. 

6.3. Composite Lya Dependences 

Through dividing our bright galaxies (R<25.5) into 
a UV-bright non-LAE and a UV-bright LAE compos- 
ite, we can look for trends based on Lya equivalent 
width. The clearest trend across these two composites 
is stronger interstellar absorption with decreasing Lya 
emission. The top panel of Figure [7] shows the highest 
Lya equivalent width quartile of S03 (top, red; WL ya = 
52.63A) and the UV-bright LAE composite (bottom, 
black; Wl yq = 50. OA) while the bottom one shows the 
second highest Lya equivalent width quartile of S03 (top, 
red; W LyQ = H.OOA) and the UV-bright non-LAE com- 
posite (bottom, black; WL ya = 8.4A). Across these 
two pairs of composite spectra, there are many interest- 
ing differences including Si II* emission, low- and high- 
ionization absorption strength, UV spectral slope, and 
velocity offset. In both the UV-bright LAE and UV- 
bright non-LAE composites, we find stronger Si II* emis- 
sion at A1265 and weaker Si II* emission at A1309 than 
the WLya = 52.63A and WL yQ = ll.OOA composites of 
S03. Furthermore, for both our composite spectra and 
that of S03, Si II* emission becomes stronger with in- 
creasing Lya emission. 

Due to poor S/N in the A > 1600A continuum, 
we only use the bluest low-ionization absorption lines 
(A1260, A1303, A1334, A1527) in determining the average 
low-ionization line equivalent widths. Figure [8] shows 
the dependences of various spectroscopic features on Lya 
equivalent width. As can be seen in the top left panel 
of Figure the average low-ionization equivalent width 
decreases from Wl/s = -1-8A to Wlis = — 1.25A for 
a change in Lya equivalent width of Wo = 8.4A to 
Wo = 50. OA. We find a similar trend in low-ionization 
absorption strength with respect to Lya equivalent width 
as S03. Note that the average low-ionization equivalent 
widths for the four S03 quartiles are different from figure 
9 of S03 as we do not include the Fe II A 1608 and the Al 
II A 1670 transitions in our measurement. For the same 
set of transitions, we find systematically stronger low- 
ionization absorption at a given Lya equivalent width. 
In the UV-bright LAE composite, the Si II A1526 tran- 
sition is only marginally detected, which differs from the 
highest Lya equivalent width quartile in S03. The de- 
gree of saturation of Si II can be measured from the ratio 
of Si II A1260 to Si II A1526, which is optically thin for 
W (1260)/Wo(1526) > 5. The ratio of Si II absorption 
line equivalent widths is W o (1260)/W o (1526) = 0.8±0.2, 
consistent with unity in the UV-bright non-LAE com- 
posite while it is W o (1260)/W o (1526) = 2.4 ± 1.5 in the 
UV-bright LAE composite. While this may indicate that 



Si II is not optically thick in the UV-bright LAE compos- 
ite, the faint nature of both transitions makes this ratio 
uncertain. With increasing Lya equivalent width, we ob- 
serve an average decrease of > 50% in the strength of the 
Si II AA 1260, 1526 and O I + Si II A1303 absorption. 
However, we note that the strength of C II A1334 instead 
increases by 20% with larger Lya equivalent width. 

Robustly determining Si IV equivalent widths is chal- 
lenging given the uncertainty in the continuum level, 
which is affected by broad stellar absorption as well as 
rcdshifted emission. In the top right panel of Figure 
Si IV AA 1393, 1402 equivalent widths are plotted for 
this work (red diamonds) and S03 (black points). The 
UV-bright LAE Si IV A1402 equivalent width measure- 
ment is likely being filled in by redshifted emission, and 
the A1393 transition is the same strength as in the UV- 
bright non-LAE composite. For this reason, we find no 
significant change in Si IV equivalent widths and con- 
clude that the transition is consistent with being opti- 
cally thin, which for Si IV corresponds to a doublet ra- 
tio of W1393/W1402 = 2. Furthermore, we do observe 
a significantly smaller C IV absorption strength in the 
UV-bright LAE composite W C iv = -1.1 ±0.35 than the 
UV-bright non-LAE composite Wciv = -2.2 ± 0.4, cor- 
responding to a decrease of — 100%. S03 reported find- 
ing 50% weaker high-ionization equivalent widths (Si IV 
and C IV) in their highest Lya equivalent width quartile 
WLya = 52.63 with respect to their other three quar- 
tiles. However, they only find a decrease of — 30% in C 
IV equivalent widths. Due to the increase in C II equiv- 
alent width with increasing Lya equivalent width, this 
may indicate that the neutral gas has a larger covering 
fraction than the ionized gas, which likely has a larger 
range of velocity offsets. Our C IV and N V profiles 
also differ qua litatively from S03 a nd high-redshift UV- 
bright SFGs (jHeckman et all 120111 ). in that we find the 
profiles to be dominated by the interstellar component, 
as opposed to the stellar component. 

As can be seen in the lower left panel of Figure [3 
for increasing Lya equivalent width, the UV contin- 
uum slopes of the composite spectra become bluer. We 
do note that the spectral slopes of the VIMOS spectra 
were scaled so that their average spectral slope matched 
that of the FORS spectra. By fitting a power law to 
the 1250 - 1600A region of the composite spectra, we 
find a UV spectral slope (Fa cx A^) of /3 = —2.25 and 
13 = —1.50 for Lya equivalent widths of Wo = 50. OA and 
Wo = 8.4A. S03 also find a bluer continuum slope for 
galaxies with stronger Lya emission. Similarly, we find 
a UV spectral slope of j3 = —1.85 and /3 = —1.40 for the 
highest and second highest Lya equivalent width quar- 
tiles from S03 respectively. The bottom left panel of Fig- 
ure[7|shows the relation between spectral slope index and 
Lya equivalent width. Although we find the same trend 
of bluer continuum slopes with increasing Lya equivalent 
width, we do find systematically bluer continuum slopes 
than S03, which may be impacted by the VIMOS flux 
calibration. Using photometrically determined /3 values, 
we find a mean spectral slope among individual galaxies 
of < p >= -1.52 ± 0.45 and < /3 >= -1.22 ± 0.47 for 
the UV-bright LAE and non-LAE composites. 

Across the UV-bright non-LAE and UV-bright LAE 
composites, we measure the average velocity offset be- 
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tween low-ionization absorption and Lya emission to re- 
main constant from Ai> em _ a b s = 640 ± 10 km s _1 to 
Aw cm _ a bs = 610 ±60 km s _1 with increasing Lya equiv- 
alent width. As is shown in the bottom right panel of 
Figure [71 our measurements are plotted as red diamonds 
while those of S03 are in black. Our results are also con- 
sistent with S03, who found velocity offset to decrease 
monotonically with increasing Lya equivalent width. We 
do note finding consistently higher velocity offsets than 
S03 at each Lya equivalent width. ISteidel et al.l (|2010f ) 
found a median velocity offset of A^em-abs = 609 km s _1 
with a scatter of 32km s _1 for their sample of 89 z^2.3 
SFGs. With the exception of baryonic and dynamical 
mass, they report not finding any significant correlations 
between velocity offset and other galaxy parameters. 

6.4. Intrinsic UV Color Dependences 

The spectral slope of the far-UV continuum is de- 
termined by the star formation rate history and the 
amount of extinction. A continuously star-forming 
galaxy has an unreddened UV spectral energy distribu- 
tion shape that remains fairly constant for ages lOMyr 
to l Gyr whose spectral slo pe, (3, ranges from -2.6 to 
-2.1 (|Leitherer et al.l 11999ft . UV-bright, 2 < z < 
3.5, st ar- for ming galaxies typica ll y have ages within 
this range (IPapovich et al.l [20011 IShaplev et ail [20011 
lErb et alil2006bl;IGawiser et al.ll2007t IKornei et al.ll2010ft 
indicating that large differences in UV spectral slope re- 
flect varying amounts of dust extintion. Due to uncer- 
tainties in the star formation histories and the form of the 
dust extinction law, we report trends based on the spec- 
tral slope index. Continuum spectral slopes are too dif- 
ficult to measure in individual spectra, however, broad- 
band photometry allows us to parameterize our data set 
based on spectral slopes as calculated from intrinsic UV 
colors. 

Figure [9] displays two composite spectra composed of 
the red (/3 > -1.4) and blue (/3 < -1.4) half of our 
UV-bright galaxies. Through fitting a power law to the 
observed photometric colors, we find spectral slopes of 
< >= -1.05 ± 0.34 and < /3 >= -1.73 ± 0.34 for 
the red and blue composites respectively where the er- 
rors represent the standard deviation. Since the spectral 
slopes of the VIMOS spectra were scaled so that their 
average spectral slope matched that of the FORS spec- 
tra, we use the photometric j3 values when examining 
the relationship among galaxy properties and spectral 
slope. Figure [TU] shows the dependence of Lya equiva- 
lent width on j3 (left panel) and the dependence of aver- 
age low-ionization equivalent width on /3 (right panel) for 
this work (red diamonds) and S03 (black points). With 
bluer spectral slope, Lya equivalent width increases from 
W = 11.0 ± 1.1 A to Wo = 27.0 ± 1.8A. Similarly, 
the average low-ionization absorption line strength de- 
creases from Wlis = -1.83A to W L is = -1.35A. We 
note that similar to the UV-bright LAE and non-LAE 
composites, while Si II absorption decreases significantly 
with bluer slopes, C II absorption does not. Addition- 
ally, the change in Lya equivalent width is not sufficient 
to fully explain this decrease in low-ionization absorption 
strength, and therefore a correlation must exist between 
spectral slope and low-ionization absorption. In contrast 
to the UV-bright LAE composite, the ratio of Si II ab- 



sorption line equivalent widths is Wo(1260)/Wo(1526) = 
1.0±0.5 and W o (1260)/W o (1526) = 0.8±0.2 in the blue 
and red composites respectively, both consistent with 
unity. A ratio of unity indicates that Si II is optically 
thick, suggesting this is the case in both composites. 

The Si IV A1393 and C IV A1549 high-ionization 
strengths do not vary substantially with spectral slope. 
With bluer spectral slope, the Si IV A1402 absorption 
line strength decreases by ~50% due to P Cygni emis- 
sion filling in the absorption. Moreover, relative to the 
red composite, the N V P Cygni profiles are more promi- 
nent. The N V A 1240 absorption equivalent width is 
W = 0.8 ± 0.3 with a large velocity offset, Av > 1200 
km s _1 . Nebular features originate from H II regions at 
the systemic redshift and indicate a larger presence of 
high mass O- and B-type stars. On the other hand, we 
do not find a C IV A 1549 P Cygni profile in the blue spec- 
tral slope composite while we do observe one in the red 
one. Finally, there is no significant change in Si II* emis- 
sion strength with spectral slope, in spite of a difference 
in Lya equivalent width. 

6.5. Redshift Evolution 

There are roughly the same number of UV-bright L AEs 
and UV-bright non-LAE objects at z~2 and z~3, al- 
lowing us to divide our sample based on redshift and 
Lya emission, generating four composite spectra. The 
low-redshift sample (2< z <2.7) is composed almost en- 
tirely of FORS objects while the high-redshift sample 
(2.7< z <3.5) is mostly VIMOS objects. In Figure 10, 
the top panel shows the two UV-bright LAE composites 
for z^2 (red, upper) and z~3 (black, lower). The bot- 
tom panel is the same as the top panel except for the 
UV-bright non-LAE composites. In this figure, the Lya 
dependences discussed in §6.31 are also apparent. These 
trends are observed at both redshifts, as the two UV- 
bright LAE composites show significantly stronger low- 
ionization absorption. The S /N of the low-redshift com- 
posites degrades quickly blueward of Lya due to poor at- 
mospheric transmission and low throughput of the spec- 
trograph below 3600 A. We present equivalent widths, 
velocity offsets and continuum spectral slopes for these 
composites as well as a z~3 all-LBG and z~2 all-BX 
composites in Table [6l 

For UV-bright LAEs, there are several interesting 
trends with decreasing redshift, including weaker Lya 
emission, stronger insterstellar absorption, stronger Si 
II* emission, and bluer spectral slopes. The change in 
Lya emission line strength is the clearest trend, as from 
z*~3 to z^2, the Lya equivalent width decreases from 
W = 66.5±6.0Ato W = 29.5±5.5A. We find interpret- 
ing this result difficult as there are a multitude of effects 
that cont ribute to the s t rength of t he Lya emission line . 
However, iNilsson et all ([20091 ) and iGuaita et ail (|2010f ) 
both find LAEs at z^2 to appear more evolved than those 
at z~3, although they do find comparable photometric 
Lya equivalent widths. Aside from Lya emission, very 
few features are seen in the z~3 UV-bright LAE com- 
posite, which may be due to the faint continua of z^3 
UV-bright LAEs resulting in a low S/N. For this reason, 
it is difficult to determine if the differences between these 
two composites reflects an evolution in galaxy properties. 
In the z~2 composite, the weak absorption strength of 
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the Si II resonant lines are again apparent as the A1260 
transition is faint and there is no A1526 absorption. In 
contrast, we do find evidence for Si II* fine-structure 
emission at AA 1265, 1533, but again we do not see Si II* 
A1309 emission. We also note that the relative strength 
of Si II* emission with respect to Si II resonant absorp- 
tion is significantly higher in our low-redshift UV-bright 
LAE composite than either the high-redshift UV-bright 
LAE composite or the high-Lya equivalent width quar- 
tile composite of S03, which was also composed of z^3 
galaxies. 

In the UV-bright non-LAE composites, we find little 
evolution in galaxy properties from z^3 to z^2 with 
the exception of stronger high-ionization absorption at 
higher redshifts. The high-ionization interstellar absorp- 
tion strengths are on average ~ 50% stronger in the high- 
redshift subsample. However, we note that the z^3 sub- 
sample is dominated by VIMOS LBGs with a larger in- 
strument FWHM. As can be seen in Tabled we find com- 
parable Lya equivalent widths, low-ionization equivalent 
widths, and velocity offsets. Therefore, this discrepancy 
may not reflect an actual evolution in high-ionization 
strength. Furthermore, there is Si II* emission at A1265 
of the z~2 composite, although we find no other evidence 
for Si II* fine-structure emission at either redshift. Thus, 
the relative strength of Si II* emission for SFGs appears 
to be stronger at z~2 than z^3. As was previously seen 
in the UV-bright LAE and UV-bright non-LAE compos- 
ites, Si II* emission is stronger with stronger Lya emis- 
sion. 

As discussed in § [31 the spectral slopes of individual 
galaxies were measured based on their photometric col- 
ors for which we found bluer spectral slopes at lower 
redshift. This trend is marginally significant among the 
UV-bright non-LAE and LAE subsets. We find bluer 
mean spectral slopes of ft-LBG-oniy — —0.97 ± 0.58 to 
fox-only = -1.39 ± 0.24 for UV-bright non-LAE galax- 
ies and /3-lbg-lae = -1.34 ± 0.33 to Pbx-lae = 
— 1.87 ± 0.51 for UV-bright LAE galaxies (errors rep- 
resent the standard deviation of the mean). The errors 
due to photometric uncertainties are comparable to the 
standard deviation of the mean. We do not attempt to 
measure the spectral slope from the continuum of the 
composite spectra as all of our high-redshift subsamples 
are composed almost entirely of VIMOS objects whose 
individual spectra were scaled to match the FORS spec- 
tra. Bluer spectral slopes at lower redshifts would in- 
dicate that UV-bright non-LAE galaxies and UV-bright 
LAEs are becoming less evolved or less dusty with cosmic 
time. While this trend of bluer spectral slopes with de- 
creasing redshift is both surprising and interesting, more 
observations will be necessary to understand if it is real. 
Finally at each redshift, we also find bluer spectral slopes 
with stronger Lya emission indicating a connection be- 
tween line and continuum extinction at both redshifts. 

7. CONCLUSIONS 

We have combined and analyzed the spectra of both 
narrowband- and continuum-selected star-forming galax- 
ies across a redshift range of 2 < z < 3.5. Through 
placing a canonical flux limit of R<25.5, we are able 
to study the relationship among Lya emission, intrinsic 
UV spectral slope, redshift, outflow kinematics, low- and 
high-ionization absorption, and nebular emission for a 



sample of galaxies spanning a similar luminosity range. 
Our results are summarized below. 

• UV-bright SFGs have spectral slope indices, /3, that 
span a large range from —2.7 to 0.1 while UV- 
bright LAEs (WLy Q > 20A) have bluer spectral 
slope indices with —2.7 < /3 < —0.8, in good agree- 
ment with model predictions for a you ng dust-free 

SFG (/3 2.6. ILeitherer et al.lll999h . Of the 59 

R<25.5 SFGs, 10 have Si II* emission and 12 have 
multiple-peaked Lya emission. Additionally, high- 
ionization absorption lines are dominated by the 
interstellar, as opposed to the stellar component. 

• UV-bright LAEs with rest-frame WL ya > 20A have 
weaker low-ionization absorption, weaker C IV ab- 
sorption, bluer spectral slopes, and stronger Si II* 
fine-structure emission than UV-bright non-LAE 
galaxies with rest-frame WL yQ < 20A. Among 
low- and high-ionization absorption strengths, in- 
dividual elements also show different dependences 
on Lya emission. 

• SFGs with bluer spectral slopes have more promi- 
nent N V AA 1238,1242 and Si IV A1402 P 
Cygni profiles, weaker low-ionization absorption, 
and stronger Lya emission. 

• From 2.7 < z < 3.5 to 2.0 < z < 2.7, SFGs 
exhibit bluer spectral slopes and stronger Si II* 
fine-structure emission. Additionally, those with 
rest-frame WL ya < 20A have comparable low- and 
high-ionization absorption strengths as well as sim- 
ilar Lya emission strengths. In contrast, there is 
a significant decrease in Lya emission strength in 
the subsample with rest-frame WL yQ > 20A, but 
we caution this may be due to selection effects. 

The observed trends with Lya emission are largely 
consistent with the physical picture put forth by S03 
where the mechanical energy from star formation and 
supernovae in high-redshift UV-bright SFGs is powering 
galaxy-scale outflows of ionized and neutral gas. Evi- 
dence for these outflows is seen as strong blueshifted ab- 
sorption as well as redshiftcd H I Lya emission relative 
to nebular features at the systemic redshift. Amongst 
the individual galaxies, a wide range of velocity offsets 
between Lya emission and interstellar absorption from 
Avcm-abs ~ 200 km s _1 to Aw cm _ a b s ~ 900 km s _1 is 
observed. In spite of this range in velocity offsets, we 
do not find a significant change in Av cm - a bs with re- 
spect to Lya emission strength in the composite spec- 
tra. We note that our velocity offsets decrease slightly 
with Lya emission and are also consistent with S03 who 
found velocity offset to decrease monotonically with in- 
creasing Lya emission strength. The velocity offset of our 
UV-bright SFG composite of Au cm _ a b s = 630 km s _1 is 
also remarkably similar to the average velo c ity off set of 
Auem-abs = 609 km s" 1 from lSteidel et all (|201d) for a 
sample of 89 continuum-selected galaxies at z ~ 2.3. 

Within the outflowing component, the high-ionization 
features have similar kinematic properties to and span 
a similar velocity range as the low-ionization features. 
However, this information says nothing abo ut the phys- 
ical distribution of the absorbing gas. iSteidel et al.l 
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(2010) find significant absorption of both low- and high- 
ionization species displaying similar radial dependences 
at galactocentric radii spanning 3 — 125kpc. This indi- 
cates that both ionized and neutral gas share not only 
the same velocity envelope, but also similar physical dis- 
tributions. 

Using the low-ionization absorption lines as a probe of 
the outflowing neutral gas, we find a direct correlation 
between Lya emission strength and the average strength 
of low-ionization absorption. This trend is most clearly 
demonstrated by Si II ions. In the outflowing component, 
a difference in line strengths indicates a change in either 
the velocity dispersion, covering fraction, or both. While 
the Si IV transition is consistent with being optically 
thin in all composites, the Si II transition and the other 
low-ionization transitions are optically thick. We do note 
that among galaxies with rest-frame WL ya > 20A, Si II 
absorption is significantly weaker and the A1526 transi- 
tion in the UV-bright LAE composite is marginally de- 
tected leading to a Si II A1260 to A1526 ratio that is 
greater than one, which wo uld be expec t ed if Si II were 
not optically thick. While ISteidel et all (|2010D found Si 
II to be unsaturated at b ~ 63kpc, those authors did find 
it to be saturated at b ~ 3kpc and b ~ 40kpc. We also 
note that due to the low S/N, the ratio of these two line 
strengths is fairly uncertain. Given the similar kinematic 
properties of different Si ionization states, different ion- 
ization state dependences on Lya emission, and similar 
physical distributions, we find it likely that the outflow- 
ing neutral gas is in the form of neutral clo uds embedded 
in ion ized gas as previously theorized by ISteidel et al.1 
(l20lQh . 

Similar to the results of S03, we find a significant cor- 
relation between bluer spectral slopes and weaker low- 
ionization absorption strength. SFGs with ages between 
lOMyr and lGyr and constant star-formation histories 
have similar UV spectral slopes, so the shape of the ob- 
served spectral slope is largely dependent on dust ex- 
tinction. UV-bright SFGs and LAEs at z~2 and z~3 are 
well-c haracterized by these ages and star formation his- 
torics dShaplev et aLll200U Eai et al.ll2008HNilsson et all 
120091 iGuaita et all 120101 ) . Therefore, we interpret this 
result as evidence for the presence of dust in the out- 
flowing component. In addition, these outflows must also 
cover a region sufficiently large to extinguish a significant 
amount of the UV continuum surface brightness. The in- 
creased prominence of nebular features for bluer /3 values 
is not surprising as bluer spectral slopes indicate either a 
lower amount of extinction, the presence of hotter more 
massive stars, or both. Each of these effects would make 
stellar photospheric features more visible. 

From a redshift of 2.7 < z < 3.5 to 2.0 < z < 2.7, 
the UV-bright non-LAE spectra look remarkably simi- 
lar. They have comparable Lya emission strengths, in- 
terstellar absorption strengths, and kinematics consistent 
with no significant evolution in SFGs with rest-frame 
WLya < 20A or the nature of their outflows across this 
redshift interval. In contrast, for SFGs with WL yQ > 20 A 
we find significantly weaker Lya emission in the lower 
redshift sample. If this trend is real, then it is consis- 
tent with UV-bright LAEs being more evolved at lower 
rcdshifts. Low- and high-ionization absorption also ap- 
pear to be stronger in the low-redshift sample, although 



due to the faint continua of z~3 UV-bright LAEs, this 
may not indicate an evolution in galaxy properties. For 
both subsets of SFGs, the lower- redshift composites have 
bluer spectral slopes, which is confirmed at marginal sig- 
nificance for the UV-bright LAE composite from the pho- 
tometric spectral slopes of individual galaxies. As there 
are many systematic effects that can contribute to the 
spectral slope, additional investigation specifically a bet- 
ter constraint on their rest-UV spectral slopes would be 
needed to conclude that lower redshift continuum- and 
narrowband-selected galaxies are less evolved, as their 
bluer colors would imply. 
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Table 1 

Subsample Criteria 



Type 


Magnitude 


w LyQ 


Nfors 


NyiMOS 


Ntotal 


UV-bright SFG 


R < 25.5 


all 


29 


30 


59 


UV-bright non-LAE 


R < 25.5 


< 20A 


15 


12 


27 


UV-bright LAE 


R < 25.5 


> 20A 


11 


18 


32 


UV-bright blue 


R < 25.5 




14 


11 


25 


UV-bright red 


R < 25.5 




11 


17 


28 


UV-faint LAE 


R > 25.5 


> 20A 


14 


8 


22 



Table 2 

AU-SFG Spectroscopic Features 



Ion 


Arcst (A) 


/ 


Wo (A) 


a (A) 


A»cm-ah S (km S !) 


LyQem 


1215.67 




19.1 


1.5 




Lyciabs 


1215.67 




-2.7 


0.3 


2400 


Si II 


1260.42 


1.007 


1.3 


0.15 


520 


O I 


1302.17 


0.0489 


2.1 


0.2 


760 


Si II 


1304.37 


0.094 


2.1 


0.2 


760 


C II 


1334.53 


0.128 


1.7 


0.2 


590 


Si II 


1526.71 


0.130 


1.2 


0.3 


580 


Fe II 


1608.45 


0.058 


0.7 


0.3 


590 


Al II 


1670.79 


1.83 


1.2 


0.4 


560 


N V 


1238.82 


0.157 


0.9 


0.2 


1320 


N V 


1242.80 


0.0782 


0.9 


0.2 


1320 


Si IV 


1393.76 


0.514 


1.5 


0.2 


550 


Si IV 


1402.77 


0.255 


0.9 


0.2 


630 


C IV 


1548.20 


0.191 


1.7 


0.3 


800 


C IV 


1550.78 


0.0952 


1.7 


0.3 


800 



Note. — Lya r and L ya^ refer to the red and blue components of the Lyc* profile. All values of Ai> are relative to Lya r . Transition 
oscillator strengths as in iPettini et all |[2002T ). W , a, and Av values listed for O I A1302 and Si II A1304 refer to the line profile for these 
two blended features assuming the rest wavelength is A = 1303. 27A. Similarly, the values for N V AA 1238.82, 1242.80 and C IV AA 1548.20, 
1550.78 refer to the blended line profiles assuming a rest wavelength of A = 1240. 8lA and A = 1549. 48A respectively. 
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Table 3 

Spectroscopic Properties 



UV-bright non-LAE UV-bright LAE UV-bright SFG 





27 


32 


59 


R 


24.46 


24.84 


24.72 


fiphot 


-1.3 


-1.5 


-1.3 


fispec 


-1.50 


-2.25 


-1.90 




640 ± 10 


610 ±60 


630 ± 20 


FWHMus 


340 ± 100 


420 ± 90 


430 ± 60 




8.4 ±0.9 


50.0 ±4.8 


19.1 ± 1.5 


w Lyab 


-4.9 ±0.9 


1.2 ±0.8 


-2.7 ±0.3 


WsiII,1193 


-1.9 ±0.3 


- 


2.0 ±0.2 


W / SiII,1260 


-1.4 ±0.25 


-1.2 ±0.2 


-1.3 ±0.15 


Woi+SiII,1303 


-2.5 ±0.3 


-1.7 ±0.3 


-2.1 ± 0.2 


Wcil,1334 


-1.5 ±0.25 


-1.8 ±0.3 


1.7 ±0.2 


WsiII,1526 


-1.8 ±0.3 


-0.5 ±0.3 


1.2 ±0.2 


WFeII,1608 


-1.1 ±0.4 




0.7 ±0.3 


WA1II,1670 


-1.5 ±0.45 


-1.0 ±0.5 


1.2 ±0.3 


M / NV,1240 


-1.0 ±0.35 


-1.1 ± 0.4 


0.9 ±0.2 


WsiIV,1393 


-1.5 ±0.3 


-1.4 ±0.35 


1.5 ±0.2 


WSUV,1402 


-1.0 ±0.3 


-0.4 ±0.35 


0.9 ±0.2 


WciV,1549 


-2.2 ±0.35 


-1.1 ± 0.4 


1.7 ±0.3 


M / OIII,1664 




0.8 ±0.45 





Note. — R and Pphot refer to the median values of the individual galaxies. /3 spec refers to the spectral slope of the composite. The 
UV-bright LAE and UV-bright non-LAE subsets contain a similar fraction of galaxies observed with FORS and VIMOS. Au cm _ a i, s and 
FWHMj/s refer to the average values of the low-ionization lines. 



Table 4 

Si II* Emission in UV Bright SFG Composite 



^rest 


A u i 


Wo 


Av 


(A) 


(10 8 s" 1 ) 


(A) 


(km s ) 


1264.74 


23.0 


0.6±0.2 


280 


1309.28 


7.00 


0.0±0.2 




1533.43 


7.40 


0.7±0.3 


140 



Note. — Properties of Si II* Emission in the UV-bright SFG composite. Eistein A— coefficients from the NIST Atomic Specta Database 
http://physics.nist.gov/cgi-bin/AtData/main_asd). The Si II* A1309 feature is not seen, so we do not report a Av. 



Table 5 

Spectroscopic Properties of z~2 and z~3 UV Bright Composites 





z~2 non-LAE 


z~3 non-LAE 


z~2 LAE 


z~3 LAE 


z~2 SFG 


z~3 SFG 


Ngal 


16 


11 


12 


20 


28 


31 


R 


24.55 


24.36 


24.69 


24.94 


24.62 


24.85 


Pphot 


-1.3 


-1.2 


-1.7 


-1.4 


-1.4 


-1.2 




640 


740 


600 


590 


630 


660 


Why aT 


8.2 ±0.7 


10.6 ±1.4 


29.5 ± 5.5 


66.5 ± 6 


12.3 ±0.8 


52.3 ±3.5 


w Lyab 


-4.1 ± 0.5 


-5.0 ± 1.5 


1.8 ±1.3 




-3.9 ±0.9 




1^3111,1193 


-2.2 ±0.4 


-1.6 ±0.6 






-2.8 ±0.4 




WsiII,1260 


-1.3 ±0.2 


-2.0 ±0.25 


-1.2 ± 0.4 




-1.3 ±0.15 


-1.2 ±0.25 


Woi+sai,i303 


-2.6 ±0.3 


-2.0 ±0.3 


-2.2 ± 0.4 


-0.8 ±0.5 


-2.3 ±0.2 


-1.1 ± 0.3 


Wcil,1334 


-1.5 ±0.25 


-1.6 ±0.25 


-2.3 ±0.5 




-1.7 ±0.2 


-1.3 ±0.3 


W / SiII,1526 


-1.7 ± 0.4 


-2.0 ±0.4 


-0.6 ±0.5 




-1.3 ±0.25 


-1.0 ±0.3 


WnV,1240 


-1.1 ± 0.5 




-1.0 ±0.3 


-0.8 ±0.4 


-1.0 ±0.25 




WsiIV,1393 


-1.4 ±0.35 


-1.5 ±0.45 


-1.6 ±0.5 




-1.7±0.25 


-1.5 ±0.35 


Wsiiv,1402 


-0.9 ±0.35 


-1.5 ±0.45 


-0.4 ±0.5 




-0.9 ±0.25 


-0.7 ±0.35 


WciV,1549 


-2.1 ± 0.4 


-3.4 ±0.45 


-1.5 ± 0.6 


-0.9 ±0.5 


-1.8 ±0.3 


-2.0 ±0.3 



Table 6 

Analogous to Table [3] We do not report equivalent widths for features not seen in the composite spectra. 
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Figure 1. Top left: the rest-frame Lyct equivalent width (Wo) distribution for our spectroscopic sample (black solid line) and our UV- 



bright LAE subsample (red dot-dashed line). All reported equivalent widths are in A. Top right: W t 



phot v - Wo for z^3 LAEs (black 

stars) and z~2 LAEs (red diamonds). Bottom left: Ai> om _ a t, s v. spectroscopic Wo where errors in Av om _ a b s are the range values found 
for different transitions. Bottom right: spectral slope index, /3 (Fx oc A* 3 ), as calculated from the photometric colors, v. Wo for R>25.5 
galaxies (cyan crosses), z~3 R<25.5 galaxies (red squares), and z ~2 R<25.5 gala xies (blue tria n gles). A typical error measurement is 
plotted in the upper right corner. LAEs from Guaita ct al. (2011) (black X's) and Nilsson et al. (2009) (black crosses) are also plotted. 
The vertical line is the traditional cut for LAEs (Wo = 20A) and the horizontal line is the division between the blue and red spectral slope 
composites. 
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Figure 2. Composite continuum-normalized rest-frame UV spectrum of 59 SFGs (R<25.5) with the composite one sigma error spectrum 
plotted in red. H I Lya emission and low- and high-ionization interstellar absorption originating from the outflowing component are the 
most prominent features. Weaker features such as stellar emission from P Cygni profiles and fine-structure emission (Si II* and C II*) are 
also visible. 
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Figure 3. A series of continuum-normalized spectra for our three subsamples: top - UV-bright non-LAE, middle - UV-bright LAE, bottom 
- UV-faint LAE. The galaxies used to construct these spectra were shifted into the rest-frame using their Lya emission line redshifts. The 
one-sigma composite error spectrum is plotted in red at the bottom of each composite. With increasing Lyo emission, interstellar absorption 
features become weaker, Si II* AA 1265, 1309 emission becomes stronger, and the spectral slopes is bluer. The A5577 skyline has been 
masked in the UV-faint LAE composite (green) at 1360A as all UV-faint LAEs lie at z = 3.1. 
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Figure 4. Zoomcd-in regions of the UV-bright LAE (top panel) and UV-bright non-LAE (bottom panel) composite spectra shown in 
Figure 3. The adjusted scale allows for a more in depth look at the interstellar absorption line profiles as well as the nebular emission 
feature Si II* AA 1265, 1309. Note the strength of the A 1265 emission as compared to the lack of A 1309 emission in both composites. The 
low-ionization transitions are all blueshiftcd with respect to Lya and the N V AA 1238, 1242 transition has a significantly higher velocity 
offset in the UV-bright non-LAE than LAE composite. The vertical lines indicate the Lya rest-frame wavelengths with respect to Lya 
emission for low-ionization absorption (dotted line), high-ionization absorption (dashed line), and Si II* emission (solid line). 
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Figure 5. Analogous to Figure 4, except showing the Si IV AA 1393, 1402, Si II A1526 and C IV AA 1548, 1550 transitions for the 
UV- bright LAE composite (top row) and UV-bright non-LAE composite (bottom row). This view shows the lack of Si IV A1402 and Si 
II A1526 absorption as well as the strength of the Si II* A 1533 emission in the UV-bright LAE composite. In the left panels, the Si IV 
redshifted emission is more apparent and the Si IV AA 1393, 1402 absorption are seen at similar velocity offsets for each composite. In the 
right panels, the decrease in Si II* A1526 emission and increase in Si II A1526 and C IV AA 1548, 1550 absorption with decreasing Lyo 
emission is clearer. 
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Figure 6. Spectroscopic incompleteness as a function of R band magnitude and spectral type is plotted for the entire catalog. Top panel 
- the black solid histogram indicates the total number of narrowband-selected objects targeted spectroscopically and the red dot dashed 
line shows the number that were spectroscopically confirmed. Bottom panel - the black solid histogram indicates the total number of 
broadband-selected objects targeted spectroscopically, the red dot dashed line shows the number that were spectroscopically confirmed, 
and the blue dashed line shows the number that were spectroscopically confirmed galaxies without Lye? in emission. Our spectroscopic 
confirmation rate is 33% for broadband-selected objects and 60% for narrowband-selected objects. The vertical dashed line represents the 
canonical R < 25.5 cut for broadband-selected galaxies. 
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Figure 7. Dependences of UV spectroscopic properties on Lya equivalent width for the UV-bright non-LAE and UV-bright LAE sub- 
samples as red diamonds, the UV-bright SFG composite in blue and the four Lyo equivalent width quartiles from S03 in black. Top left - 
average low-ionization equivalent widths (Si II A 1260, O I + Si II AA 1302,1304, C II A 1334, Si II A 1526). Top right - Si IV equivalent 
widths. Bottom left - spectral slope, /3, as measured from the composite spectra. The error bars for our /3 values reflect the error in 
spectral slope index and do not take into account the uncertainty in the artificial flux calibration of the VIMOS spectra. Bottom right - 
velocity offset Au em _ abs . The errors in velocity offset represent the range of velocity offsets for different low-ionization lines. The trends 
in spectroscopic properties with increasing Lyo emission originally observed by S03 are confirmed for our smaller sample of SFGs spanning 
a larger range in redshift. 
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Figure 8. Top - the highest Lya equivalent width quartile composite (WL ya = 52.63A) from S03 in red and the UV-bright LAE composite 
(Wi,y a = 50. OA) in black. Bottom - the second highest Lya equivalent width quartile composite (WLya = ll.OOA) from S03 in red and 
the UV-bright non-LAE composite (W-^ ya = 8.4A) in black. Our composite spectra show similar trends in both interstellar absorption 
and Si II* emission as those of S03. In the top panel, S03 find stronger Si II* A1309 emission, weaker Si IV A1402 nebular emission, and a 
more prominent C IV AA 1548, 1550 profile. Note that the average low-ionization equivalent widths for the four S03 quartiles are different 
from figure 9 of S03 as we do not include the Fe II A1608 and the Al II A 1670 transitions in our measurement. 
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Figure 9. Composite spectra of SFGs divided into two subsamples based on intrinsic UV spectral slopes as calculated from photometric 
colors. Top - the red UV continuum slope composite (fiphot > —1.4). Bottom - the blue continuum slope composite ((3 p hot < —1.4). 
With bluer spectral slope, Lyo emission increases, low-ionization absorption become weaker and there is significantly stronger Si IV A1402 
nebular emission. 



30 



20 



3 <n 



10 



in 







3 



2 



-1 
P 







2 







Figure 10. Left - the dependence of Lyo equivalent width on spectral slope, 0, when the dataset is subdivided based on spectral slope 
for this work (red diamonds) and S03 (black points). Right - the dependence of low-ionization absorption on spectral slope. The trends in 
decreasing Lyo emission and increasing low-ionization absorption for bluer spectral slopes observed by S03 are confirmed. 
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Figure 11. Composite spectra of UV-bright LAEs and UV-bright non-LAE objects for 2.0 < z < 2.7 BXs in red and 2.7 < z < 3.5 LBGs 
in black. Top - a composite of 13 BX-LAEs is plotted in red above a composite of 21 LBG-LAEs in black. For our sample of galaxies, 
few features are seen in the z~3 composite, however we do see a decrease in Lya emission with redshift. In spite of the low S/N of the 
BX-LAE composite, the relative strength of Si II* AA 1265, 1533 emission to resonant Si II absorption is significantly higher than any other 
composite. Bottom - the BX-only composite composed of 15 BXs is plotted in red above the LBG-only composite consisting of 11 LBGs in 
black. These composite spectra are quite similar and show no significant evolution of spectroscopic properties across this redshift interval. 



